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INTRODUCTION 



A . BACKGROUND 

Reductions in size, weight and cost have become 
increasingly important design considerations in the 
development of every component destined for use aboard United 
States Navy vessels. Significant reductions in weight have 
already been achieved for main propulsion machinery with the 
advent of the gas turbine propulsion plant. Weight reduction 
is also possible through redesign of auxiliary systems and 
other large components such as main steam condensers. The use 
of heat transfer enhancement techniques for condenser tubing 
has been one avenue receiving considerable attention in the 
pursuit of this goal. Smaller more efficient condensers would 
result in reductions in capital cost, space requirements, 
weight and operating cost. 

The four primary factors that affect the performance of 
condensers are the thermal resistances on the coolant side, 
across the tube wall, the vapor side, and also that due to 
fouling. Each can be reduced to minimize its effect, which in 
turn enhances the overall heat transfer rate of the system. 
In experimental studies, the fouling resistance can be 
neglected if the tube is thoroughly clean prior to testing. 
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Enhancement techniques can be categorized into two major 
groups, active and passive. Guido [Ref. 1] compiled a fairly 
comprehensive listing of both groups of enhancement techniques 
for the tube side heat transfer rate. Passive techniques such 
as internal helical ribbing and displaced promoters (the use 
of wire mesh or HEATEX inserts) will be discussed in detail in 
a later section. It should be noted, however, that both of 
these techniques have a corresponding increase in pumping 
power requirements. Wall thermal resistance can be reduced 
passively by decreasing the wall thickness or changing the 
tube material. Similar to enhancement techniques for the tube 
side, the shell side resistance can also be reduced by both 
active and passive means. The passive enhancement techniques 
of wire- wrapping and tube corrugation will also be covered in 
a later section. 

A variety of corrugated or roped tubes are available for 
use in condenser applications. Figure 1 shows the profile of 
one type of corrugated tube which has a groove radius smaller 
than the ridge radius. Wolverine Tube Company manufactures a 
low pressure drop (LPD) and a maximum heat transfer (MHT) tube 
of this type (trade name: Korodense) , available for commercial 
use. The MHT tubes have a larger groove depth compared to the 
LPD tubes [Ref. 2] . Other types of corrugated tubes are also 
available which have a larger groove radius than the ridge 
radius. Marto et al . [Ref. 3] present test results on various 
other corrugation geometries. 
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B. CONDENSATION APPLICATIONS FOR ENHANCED TUBES 

In 1971 Withers and Young [Ref. 4] performed a theoretical 
analysis on a power plant main condenser. Their findings 
indicated a substantial saving in almost every condenser 
aspect if the smooth tubes were replaced by Korodense roped 
tubes. Similar results were obtained by Webb et al . [Ref. 5] . 
They indicated a possible 30 to 40 percent increase in the 
overall heat transfer coefficient when analyzing a tube-for- 
tube replacement of smooth to Korodense tubes. This benefit 
was calculated for two separate power plants of different heat 
duty. An ongoing comprehensive study of 14 Tennessee Valley 
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Authority (TVA) power plant condensers, retubed with 
corrugated tubes, has yielded significant economic and 
performance benefits [Refs. 6 and 7] . Of the 14 condensers 
retubed by TVA, their Gallatin Unit 1 has provided long-term 
results with over 12 years of successful operation since 
retubing. In 1983 Boyd et al . [Ref. 8] studied the 
performance of TVA's Gallatin Unit 1 plant after only three 
years of operation. Their findings yielded up to a 43% 
increase in the overall heat transfer coefficient when clean 
tubes were analyzed. However, data also indicated that the 
corrugated tubes fouled quicker than plain tubes. Rabas et 
al. [Ref. 9] studied 12 of the retubed TVA condensers and 
found fouling rates for the enhanced tubes ranging from about 
the same as to almost twice that experienced in plain tubes. 
Although the fouling rates were higher, Rabas et al. [Ref. 9] 
reported the fouling resistance values for the LPD Korodense 
tubes never exceeded the minimum allowable TEMA standard. 

Union Carbide [Ref. 10] performed analysis on roped 
condenser tubes with a Thin Film Condensing Promoter 
(T.F.C.P.), which is essentially a wire-wrapped roped tube. 
Their findings indicated an increase in the heat transfer 
coefficient by a factor of between 2 to 3 times that found 
with a smooth tube under condensate loaded conditions. 
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C. NAVAL POSTGRADUATE SCHOOL CONDENSATION RESEARCH 

Titanium has a number of properties which make it a 
desirable alternative to the conventional copper- nickel tubing 
currently in use today on board Naval vessels. Its high 
strength to weight ratio enables the tubes to be designed with 
a much thinner wall thickness, while still maintaining the 
durability of the overall system. Titanium also has a high 
resistance to erosion and sea water corrosion. This leads to 
reduced shutdown time for tube cleaning and replacement. One 
of the primary disadvantages of titanium tubing is its poor 
thermal conductivity. Titanium has approximately 50% the 
thermal conductivity of copper-nickel [Ref. 11]. This 
disadvantage can be partially offset by a decreased wall 
thickness and also through the use of enhanced surfaces. The 
Naval Postgraduate School (NPS) , in conjunction with the Naval 
Surface Warfare Center, Carderock Division, Annapolis 
Detachment, has been testing the enhancement characteristics 
of various types of titanium tubing. 

In order to examine the heat transfer characteristics of 
various types of enhanced tubes, a suitable test platform must 
be utilized. While the ideal test platform would be a full 
scale model of the actual condenser used in the field, the 
cost associated with such an apparatus would be prohibitive. 
The condensation test platform, constructed by Krohn [Ref. 12] 
at NPS, has been in use since 19 82. Although several 
modifications have been made by various researchers since that 
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time, the overall apparatus and the operation remain basically 
unchanged. Van Petten [Ref. 13] lists a comprehensive summary 
of the work completed by eight NPS researchers from 1982 
through 1987. Most of the work concentrated on the 
performance of various types of externally finned tubes. 
Mitrou [Ref. 14] , in addition to studying finned tubes, also 
performed experiments on wire-wrapped smooth copper tubes. 
More recent work conducted by Coumes [Ref. 15], Guttendorf 
[Ref. 16] and Swensen [Ref. 17] dealt with further studies of 
finned copper tubes. O'Keefe [Ref. 18] studied the 
performance of various titanium tubes. His tests included 
smooth tubes, wire- wrapped smooth tubes, LPD Korodense tubes, 
and wire-wrapped LPD Korodense tubes. Other researchers 
including Kanakis, Brower and Mazzone [Refs. 19-21] have 
studied tube bundle configurations for smooth and corrugated 
tubes, both with and without wire-wrapping . 

D . OBJECTIVES 

The main objectives of this study are as follows: 

1. To verify the repeatability of the calculated 
coefficients for both coolant side and vapor side heat 
transfer correlations obtained by previous researchers on 
wire-wrapped smooth and corrugated titanium tubes. 

2. To manufacture and collect condensation data on smooth 
and MHT Korodense titanium tubes with varying fractional 
coverage due to wire -wrapping. 

3. To determine if any relationship exists for optimal 
fractional coverage or pitch to diameter ratio between 
wire-wrapped smooth tubes and wire-wrapped MHT and LPD 
Korodense tubes . 
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II. LITERATURE SURVEY 



A. INTRODUCTION 

Condensation on horizontal tubes has been an intensely 
studied heat transfer mechanism throughout the years. New 
developments in condenser tubing require continued 
experimental and theoretical research in this field. 

There are two fundamental modes of condensation on 
horizontally mounted condenser tubes, filmwise and dropwise. 
Dropwise condensation can produce heat transfer rates as much 
as an order of magnitude larger than those associated with 
filmwise condensation. Various surface coatings have been 
developed in order to promote dropwise condensation and take 
advantage of the higher heat transfer rates. However, with 
time these coatings lose their effectiveness and filmwise 
condensation eventually returns. For this reason, condenser 
design calculations are normally based on the more 
conservative filmwise mode of condensation [Ref. 22] . 

Internally, the primary factors affecting thermal 
resistance are the intensity of turbulence and the thermal 
conductivity of the coolant medium. The degree of turbulence 
can be changed by altering the internal tube wall geometry or 
through the use of displacement promoters. 
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B. VAPOR SIDE CONSIDERATIONS 



Condensation on horizontal tubes consists of a complex two 
phase process. A suitable theoretical model must be utilized 
in order to accurately predict these phenomena. 

1. Horizontal Smooth Tubes 

Nusselt [Ref. 23] developed the foundation for the 
study of filmwise condensation on horizontal tubes in 1916. 
His work was formulated for a quiescent vapor condensing on a 
single horizontal tube. Due to the increase in the thickness 
of condensate as gravity forces it around the sides of the 
tube, the local heat transfer coefficient decreases around the 
tube circumference. Nusselt developed a theory for the 
average heat transfer coefficient around the tube, accounting 
for the lower resistance at the top of the tube where the film 
was the thinnest and the higher resistance at the bottom of 
the tube where the film was the thickest. As a result, he 
developed the following equation: 



h 



O 



0.728 



k f 3 9 P t (Pf-Py) 
V-f D o (Tsat-Tn) 



1/4 



( 2 . 1 ) 



where the fluid properties are evaluated at the film 
temperature, given by: 



Tf = 





( 2 . 2 ) 



Nusselt 's theory has been extensively studied since 
1916, and with the imposed assumptions, has been found to be 
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generally valid [Refs. 24, 25] . It has also been found to be 
accurate for cases which do not conform to Nusselt's original 
assumptions, such as variable wall temperature [Ref. 26] . One 
of the major problems encountered in applying Nusselt's theory 
in the design of condensers arises from his assumption of a 
quiescent vapor. While in theory, and in some limited 
practical applications, the assumption of a stationary vapor 
can be justified, most steam condensers operate under 
conditions where the vapor is travelling at some constant 
velocity . 

With the vapor in motion, sufficient shear forces can 
be developed which significantly thin the condensate film and 
increase the outside heat transfer coefficient. Early 
theoretical work by Shekriladze and Gomelauri [Ref. 27] , which 
accounted for vapor shear, resulted in the following equation 
to approximate the mean Nusselt number. 



— = 0.64(1 +(1 + 1 . 69F) 1/2 ) 1/2 



( 2 . 3 ) 



where : 



F = 



9 D 0 \i f h 



o_rf 11 fg 



( 2 . 4 ) 




and the two phase Reynolds number is given by: 




( 2 . 5 ) 
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Lee and Rose [Ref. 28] compared vapor shear models with 
experimental results and found that the Shekriladze-Gomelauri 
results were more conservative than more rigorous developments 
performed by other researchers due to their simplified 
approximation for the interfacial shear stress. Fujii et al . 
[Ref. 29], in a more recent study, performed an extensive 
development on the outside of a horizontal tube. From 
experimental data, they also developed an empirical 
formulation for the condensation of steam on a horizontal tube 
which includes vapor velocity effects: 



Nu 



Re 1/2 



= 0.96 F 1/5 



( 2 . 6 ) 



where F and Re 2( p are defined in equations (2. 4) and (2.5) . For 
situations where the surface shear forces dominate, equation 
(2.6) tends to more accurately predict the vapor side heat 
transfer coefficient for steam. 

As discussed previously, the primary modes of external 
tube enhancement covered in this study are wire- wrapping and 
the addition of helical grooves. Both methods have associated 
advantages and disadvantages . 

2. Wire-Wrapped Smooth Tubes 

The benefits of wire- wrapping smooth tubes have been 
attributed to two primary effects: the thinning of the 
condensate film between the wires due to surface tension 
forces and the better drainage of the condensate along the 
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wire, (Fujii et al . [Ref. 30] and Marto et al . [Ref. 31] ) . 
Thomas [Ref. 32] conducted tests on vertical tubes with wires 
attached vertically along the tubes and demonstrated this 
latter enhancement effect. Horizontal smooth tubes with 
coiled wire were later studied by Thomas et al . [Ref. 33] and 
their findings indicated a two fold increase in the heat 
transfer coefficient over that predicted by Nusselt theory for 
a smooth tube. Marto et al . [Ref. 31] determined an optimum 
wire pitch to diameter ratio between 5 and 7 for steam 
condensation. They achieved an enhancement between 1.5 and 
1.8 for the heat transfer coefficient over plain tubes. 

For low surface tension fluids (R-ll and ethanol), 
Fujii et al . [Ref. 30] obtained results 2 to 3 times those 
predicted by Nusselt theory, with an optimum wire pitch to 
diameter ratio of 2. O'Keefe [Ref. 18] compared the 
enhancement ratio versus fractional coverage (the percentage 
of the tube's surface covered by wire) on a titanium smooth 
tube. Together with data obtained by Mitrou [Ref. 14] and 
Sethumadhavan and Rao [Ref. 34], he determined that the 
optimum fractional coverage was between 10 and 30%. This 
fairly broad range of values was attributed to the possible 
effects of the different tube material thermal conductivities 
used in the studies . 
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3 . Corrugated Tubes 

In 1971, Withers and Young [Ref. 35] performed a 
design study of two condenser bundles with heat duty of 25 
million Btu/hr, one with bare (smooth) tubes and the other 
with corrugated (roped) tubes. Coolant flow pressure drop was 
one of the matching parameters for the comparison. This 
requirement made it necessary to use a larger corrugated tube 
(1.00 inch outside diameter (OD) ) , than the smooth tube (0.75 
inch OD) . They found an enhancement of 88% for the outside 
bundle heat transfer coefficient. Catchpole and Drew (Ref. 
36] tested single corrugated tubes and determined that the 
steam side heat transfer coefficient increased with increasing 
groove density and passed through a maximum with increasing 
groove depth. Their conclusion was that the smaller the helix 
angle (measured from the tube axis) , the worse the condensate 
drainage, resulting in a reduction of the outside heat 
transfer coefficient. Similar results were reported by Marto 
et al . [Ref. 3] and Mehta and Rao [Ref. 37] . Eleven different 
corrugated tubes of three primary types were tested by Marto 
et al . [Ref. 3] . Within each corrugation type, the largest 
value for the outside heat transfer coefficient was obtained 
at small pitch to diameter ratios (P/D) . As P/D increased, 
the enhancement ratio decreased. In tests conducted by Mehta 
and Rao [Ref. 37], they found the condensing coefficient 
enhancement ratio decreased from 1.2 6 to 1.1 when the 
corrugation groove pitch was increased from 6.35 mm to 12.5 mm 
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with a constant groove depth of 0.63 mm. They also tested the 
effect of holding the groove pitch to groove depth ratio 
constant. A maximum enhancement ratio of 1.38 was achieved 
with the smallest groove depth and pitch. Zimparov et al . 
[Ref. 38] confirmed the same trend, using 27 mm outside 
diameter corrugated tubes, with the enhancement ratio 
increasing from 0.99 to 1.22 for a corrugation groove pitch 
decreasing from 16.9 mm to 6.5 mm. 

Cunningham and Bahjernejad [Ref. 39] examined the 
effects of condensate inundation on corrugated tube bundles. 
They tested two different multistart tubes, one with two 
helical starts and the other with six starts. Their results 
indicated severe limitations in performance due to inundation 
effects for tube bundles more than a few rows deep. 
Additional inundation tests were conducted on multistart tubes 
by Cunningham and Boudinar [Ref. 40] . They reported similar 
results as those of Cunningham and Bahjernejad [Ref. 39] , and 
additionally found that the more severe the indentation 
(groove depth) on the roped tube, the greater was the 
reduction in vapor side heat transfer for the tube bundle. 

Rabas [Ref. 41] studied an extensive data bank of 
results published by numerous researchers. He found that the 
enhancement in the outside heat transfer coefficient was 
negligible for tubes with grooved geometries that are suitable 
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for ball cleaning 1 of the inside surface. He also compared 
Che collected data with the results predicted by three 
different correlations. The Wolverine correlation [Ref. 2]: 



h Q = 0.725 C n 



k f 2 p/ g h 



1/4 






n t Pf D o < T sac~ T wo ) 



where : 



C n = Z* [n t ] w 



( 2 . 7 ) 



( 2 . 8 ) 



Z* 



38 (o ) 0-65 + 



( 2 . 9 ) 



and: 



PW = [(5 .31E-7) + [0.195 (o>) ] 1/6 



( 2 . 10 ) 



• Values for a> and 9 * are listed for specific tube 
geometries by Wolverine in Reference 2. 

• The variable C n in Equation (2.7) must be applied using 
English units . 



and the correlation developed by Mehta and Rao [Ref. 37] : 



(kf 2 P t g/v-t) 



2\l/3 



= 1.54 e 



2058 



4r 



-O^e " 7 - 519 



( 2 . 11 ) 



for 0 < 6 < 15xl0" 4 , and: 



■'•ball cleaning is a method of removing deposits from the 
inside of condenser tubes by periodically recirculating sponge 
rubber balls through the tubes with the coolant. 
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= 1.92 e' 13 - 20 



ill 



-0 . 33e ' 



( 2 . 12 ) 



(k f 2 p/ a/ H f 2 ) 1/3 

for 15xl0' 4 < 0 < 200xl0' 4 , were both found to overpredict the 
outside heat transfer coefficient. Values similar to those 
obtained from the experimental data were predicted by the 
Yorkshire correlation [Ref. 42] given by: 



k 0 = c 8 



V P t 2 9 

Pr r 



1/3 



( 2 . 13 ) 



where : 



Cg 1 






[ L °l 


p' 


D o\ 


Do. 



-0.75 



( 2 . 14 ) 



and 4> c varies with tube profile and can be found in Reference 
42. However, Rabas [Ref. 41] cautioned its use due to the 
results being close to unity. In conclusion, he recommended 
using smooth tube correlations and neglecting the minimal 
enhancement associated with corrugated tubes for design 
calculations . 

Dreitser et al . [Ref. 43] determined a correlation for 
the enhancement ratio of corrugated tubes based on three 
primary geometric ratios: 1) The ratio of the radius of 
curvature of the corrugation ridge to the outside diameter, 
R/D 0 ; 2) the ratio of groove pitch to the outside diameter, 
P/D 0 ; and 3) the ratio of the external groove diameter to the 
outside diameter, d 0 /D 0 . 
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exp 


3.65 1— -2 
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D o . 



(2.15) 



They found chat if any of the three ratios were increased 
while keeping the others constant, the enhancement ratio would 
decrease . 

Cunningham and Fidler [Ref. 44] made a comparative 
study of the effects of vapor shear on multistart corrugated 
tubes and smooth tubes . They found that the corrugated tubes 
maintained their enhancement advantage over the smooth tube 
when subjected to increased vapor velocity. 

4. Wire -Wrapped Corrugated Tubes 

Wire -wrapping of corrugated tubes has been shown by 
various researchers to significantly reduce the effect of 
condensate inundation (Memory et al . [Ref. 45] , Union Carbide 
[Ref. 10] , and Marto and Wanniarachchi [Ref. 46] ) . 

Union Carbide [Ref. 10] tested steam condensation for 
a wire-wrapped roped tube in a bundle arrangement. They 
reported nearly identical performance for both the unloaded 
and the condensate loaded conditions tested. Condensing 
coefficients approximately 2.5 times that of plain tube 
results were achieved for the condensate loaded condition. 
Memory et al . [Ref. 45] obtained a 60 to 68% enhancement over 
smooth tube results when condensing R-113 on a wire- wrapped 
roped tube bundle arrangement containing four tubes in a 
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vertical column. Their tests indicated an optimal wire pitch 
to diameter ratio of 7.5. 

O'Keefe [Ref. 18] tested the enhancement effects of 
three different wire diameters (1.6 mm, 1.0 mm, and 0.5 mm) 
when used to wrap an LPD Korodense tube. Enhancements ranging 
between 20 and 30% above smooth tube results were obtained 
with the maximum obtained for the 0.5 mm wire with a wire 
pitch to wire diameter ratio of 14. 

C. COOLANT SIDE CONSIDERATIONS 
1. Smooth Tubes 

Numerous correlations have been developed for use in 
the analysis of the coolant side heat transfer correlation 
associated with the turbulent flow regime inside a pipe (Re > 
10,000) . Many of the earlier correlations have been presented 
in the form: 

Nu = C i Re m Pr n ( 2 . 16 ) 

This form has been used for several well-known 
correlations including those developed by Dittus and Boelter 
[Ref. 47] : 

Nu = 0.023 Re°' e Pr°- i ( 2 . 17 ) 

and Colburn [Ref. 48] 
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Nu = 0.023 Re°- a Pr 1/ 3 



( 2 . 18 ) 



A correction factor for equation (2.18) was developed 
by Sieder and Tate [Ref. 49] : 

Nu = 0.027 Re°‘ e Pr 1/2 ( 2 . 19 ) 

This modification compensates for the variation in the coolant 
viscosity when large temperature differences exist between the 
bulk coolant and the inner tube wall temperature. All fluid 
properties for equations (2.17), (2.18), and (2.19) are 

evaluated at the mean coolant bulk temperature, T m . 

T = (Tl * Til ( 2 . 20 ) 

m 2 

In more recent years, Sleicher and Rouse [Ref. 50] and 
Petukhov [Ref. 51] developed equations (2.21) and (2.24) 
respectively, which are applicable over a wider range of 
Prandtl numbers. 

Nu = 5 + 0.015 Re ( c Pr w d ( 2 . 21 ) 

where : 

c = 0.88 - 0,24 ( 2 . 22 ) 

4 + Pr„ 



and, 



d = + 0.5 exp (-0.6 Pr w ) 



( 2 . 23 ) 



Nu 



(e/8) Re Pr 

K t + it, (e/8) 1/2 (Pr 2/3 -l) 



( 2 . 24 ) 
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where : 



e = (1.82 log (Re) - 1.G4)' 2 (2.25) 

/q = 1 + 3 .4 e (2.26) 

K 2 = 11.7 +1.8 Pr' 1/2 (2.27) 

Lorenz et al . [Ref. 52] compared their experimental results 
for cold water (i.e. large Prandtl number) with seven of the 
more common inside heat transfer correlations. They found 
that the Petukhov- Popov and Sleicher- Rouse correlations agreed 
to within 5% of the experimental data. The other correlations 
typically underpredicted the data by between 5 and 15%. 

2 . Corrugated Tubes 

The internal ridges of the tube wall produced when 
manufacturing corrugated tubes cause two primary effects. The 
irregularity of the wall surface results in an increased 
frictional effect and causes the flow to swirl, thereby 
requiring additional pumping power to achieve flow rates equal 
to similar diameter smooth tube flow. The increase in 
turbulence also has the effect of increasing the inside heat 
transfer rate. In the study by Withers and Young [Ref. 35], 
the total inside enhancement ratio of the tube bundle tested 
was 56% above that obtained for a smooth tube bundle designed 
with an equal pressure drop. Catchpole and Drew [Ref. 36] 
found the inside heat transfer coefficient to increase with 
increasing groove depth and groove frequency. However, they 
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stressed that the increase in performance was obtained at the 
expense of significantly greater friction factors and therefor 
increased pumping power and operating cost. In comparing 
their results with the Prandtl analogy [Ref. 53] they found 
that, due to the turbulence caused by the presence of grooves, 
the inside heat transfer coefficient followed the trend 
predicted by the Prandtl analogy. 

Mehta and Rao [Ref. 37] presented their findings in 
terms of the severity factor, the ratio of groove depth 
squared to the product of groove pitch and inside diameter. 
This dimensionless group was chosen because it has been found 
to be a useful criterion for evaluation of heat transfer and 
frictional characteristics of plain and spirally enhanced 
tubes. An improvement of 260% was obtained above smooth tube 
results with the higher inside heat transfer coefficient 
corresponding to higher severity factors. Of the eleven tubes 
tested by Marto et al . [Ref. 3], inside enhancement was 
reported for all the tube geometries examined, with the 
magnitude ranging from 1.6 to 5.1 times a smooth tube. 
Zimparov et al . [Ref. 38] tested nine tubes with differing 
corrugation geometries and obtained results similar to 
previous researchers. They achieved inside enhancement ratios 
between 1.77 and 2.73, with an increase in the friction factor 
of 2 to 4 fold. 

For corrugated tubes, Withers [Ref. 54] modified 
existing correlations and developed the following: 
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Nu = 



y/778 Re Pr 

7 . 22(P/D i )' l/l /PF [(e / D i )Re y /I78\ 



0 . 127 



( 2 . 28 ) 



where : 
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+ 3.75 






v > 





( 2 . 29 ) 



and: 




N 8 



1 

2.46 In [r + (7 /Re) *] 



( 2 . 30 ) 



The variables m and r were determined for each tube by a 
curve- fitting procedure. He obtained agreement between 
experimental and predicted data to within 10%. Wolverine 
currently utilizes a form of this correlation in their 
Engineering Data Book II [Ref. 2] . 

In a comparison of the LPD and MHT Korodense tubes 
manufactured by Wolverine [Ref. 2], they report a 25% higher 
inside heat transfer coefficient for the MHT Korodense tube 
than for the LPD Korodense tube. The MHT Korodense tube also 
has twice the pressure drop of the LPD Korodense tube. 

After compiling an extensive collection of published 
data, Rabas et al . [Ref. 55] recommended a minor modification 
to the leading coefficient in the Yorkshire correlation [Ref. 
42] to better fit the experimental data. The Yorkshire 
correlation is similar in form to equation (2.17) with an 
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additional complex leading coefficient (E 1 ) to compensate for 
the corrugation tube geometry. 

Nu = 0.027 E 1 Re 0 - 8 Pr 0 4 ( 2 . 31 ) 

where : 



= 1.0 + 1.182 E 2 Ej (P/D^ 0 - 1 


( 2 . 32 ) 


E 2 = cos ( a ) 


( 2 . 33 ) 


a = tan ~ 1 (P/nD i ) 


( 2 . 34 ) 


ln(E 2 ) = - t a j [ln(e/P i ) F' 
j~o 


( 2 . 35 ) 



and : 

a 0 = 5.0766, a 1 = 4.4224, a 2 = 1.1240, a i = 6.700X10' 2 

The modification by Rabas et al . [Ref. 55] changed the 
exponent in equation (2.32) from 0.7 to 0.406. 

Panchal and France [Ref. 56] conducted a study of 
large pitch spirally indented tubes for possible application 
to pressure drop sensitive systems. Their data resulted in 
internal enhancements as high as 100% over smooth tube data 
and they reported agreement within a 20% error band for the 
results predicted by the correlation proposed by Mehta and Rao 
[Ref. 57], given by: 
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Nu 



(2.36) 



v/772 Re Pr 

P Pr 0,55 (Re\flj2) 0,11 + y 

whe re : 

P = 9 . 4 (P/D,) '°' 15 (e/Dj 0,11 (2.37) 

and y is defined as in equation (2.29) . It should be noted 
that in equation (2.37) is constant for any given tube 
configuration . 

3 . Tubes With Inserts 

Numerous types of inserts are available for use in 
enhancing the inside heat transfer coefficient. Some of the 
types used by previous researchers at NPS include: twisted 

tape, wire-wrap, and HEATEX (a wire mesh insert supplied by 
CAL GAVIN) . A number of correlations have been developed for 
determining the inside heat transfer coefficient when using 
these various inserts [Refs. 58-60] . These correlations for 
HEATEX, however, are for specific loop densities of the wire 
matrix and are not universally valid. HEATEX inserts not only 
significantly increase the overall heat transfer rates, but 
also inhibit fouling. Gough and Rogers [Ref. 61] presented 
results collected from a condenser fitted with HEATEX inserts 
which was in operation for over four years without being shut 
down for cleaning. Their findings indicated the actual 
fouling resistance was over 96% less than the allowable limits 
after a service life of 26 times that of the original unit 
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with no insert. However, it should be noted that the use of 
an insert has a pronounced effect on increasing the pressure 
drop across the tube. Inserts are used in experimental 
applications to realize a larger temperature drop in the 
coolant flow, thus decreasing the uncertainty of the heat 
transfer calculations. 
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III. EXPERIMENTAL APPARATUS 



A. SYSTEM OVERVIEW 

The apparatus used for this research was originally- 
constructed by Krohn [Ref. 12] . The initial design has been 
modified by several previous researchers (in particular 
Swensen [Ref. 17]) and has been used successfully by Coumes 
[Ref. 15], Guttendorf [Ref. 16], Swensen [Ref. 17] and O'Keefe 
[Ref. 18] to test condensation on single horizontal tubes of 
various configurations. Figure 2 contains a general schematic 
of the overall system. The boiler, a 0.30 meter diameter 
Pyrex glass cylinder, contains ten 4 kW, 440 Volt Watlow 
immersion heaters which are used to generate steam from 
distilled water. From the boiler, steam passes up through a 
cylindrical section of Pyrex glass with an inside diameter of 
0.15 meters and a length of 2.13 meters. Two 90 degree Pyrex 
glass elbows redirect the steam back down a second similar 
cylindrical section of 1.52 meters in length. The steam then 
enters the stainless steel test section containing the 
horizontal condenser tube mounted as shown in Figure 3. The 
test section also contains a circular glass viewing port so 
the condensation process can be observed while conducting 
tests. Excess steam, not condensed on the sample tube, passes 
through the test section and into the final Pyrex glass 
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Figure 2. Schematic of the Single Tube Test Apparatus 
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Figure 3 




Schematic of the Test Section 
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cylinder containing an auxiliary condenser. The auxiliary 
condenser section collects all the condensate, where it 
returns to the boiler via a gravity drain in the baseplate. 

Cooling water for the single horizontal tube in the test 
section was supplied by two centrifugal pumps connected in 
series and drawing suction from an adjacent cooling water 
sump. A throttling valve and calibrated flowmeter allowed for 
control of the cooling water flow rate. The auxiliary 
condenser cooling water was supplied directly from the 
building water main. After setting the desired coolant flow 
rate through the single horizontal tube, the saturation 
temperature in the test section was controlled by adjusting 
the throttle valve on the auxiliary condenser coolant flow 
line . 

A suction port at the base of the auxiliary condenser 
section allowed for the removal of non- condensible gases. 
Figure 4 shows a schematic of the associated system. Non- 
condensibles drawn from the auxiliary condenser passed through 
an internal condensing coil where, any associated steam that 
might be present was condensed and collected in a plexiglass 
vacuum chamber for later removal. The non-condensibles then 
passed through the vacuum pump and were expelled to the 
atmosphere . 
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Figure 4. Schematic of the Purging System and Cooling Water 
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ORAIN 



B . SYSTEM INSTRUMENTATION 



The 440 Volt immersion heater power supply was controlled 
by a silicon controlled rectifier mounted on the electrical 
switchboard. Poole [Ref. 62] details the power calculation 
for the input into the data acquisition system. 

Three separate methods were used to monitor the pressure 
of the system: 

1. Setra model 204 pressure transducer. 

2. System saturation temperature converted into pressure. 

3. Heise solid front pressure gage (visual reference only) . 

Temperature measurements were obtained using several 
devices. The vapor temperature of the system was measured by 
a Teflon coated and a metal sheathed type-T copper/constantan 
thermocouple located just upstream of the horizontal tube 
being tested. Condensate return and ambient temperatures were 
also measured with Teflon coated type-T copper/constantan 
thermocouples. The inlet and outlet coolant temperature of 
the condenser test tube were measured with three separate 
probes : 

1. A Teflon coated type-T copper/constantan thermocouple. 

2. A ten junction, Teflon coated type-T copper/constantan 
thermopile . 

3 . An HP 2804A quartz crystal thermometer. 
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The three temperature probes at the outlet of the horizontal 
tube were located downstream of a mixing chamber in order to 
obtain a representative average coolant discharge temperature. 
Figure 3 contains the locations of these temperature probes. 

An HP-3497A data acquisition system, linked to a HP-9826A 
computer, was used to collect the voltage readings from all 
the measurement devices. The data collection and reduction 
program "DRPML", listed in Appendix C, allowed for the 
processing and storage of all the data by the HP-9826A 
computer . 

Swensen [Ref. 17] provides calibration data for the 
pressure transducer, thermocouples, quartz thermometer and 
flowmeter . 

C. TUBES TESTED 

Eight enhanced tubes were manufactured for this thesis. 
Tests were also conducted on tubes studied by previous 
researchers to ensure repeatability of results. Table I 
contains a listing of all the tubes tested and additional 
tubes used in the comparison of results. All of the tubes 
assembled for this study were made from a single titanium MHT 
Korodense tube manufactured by Wolverine. Seven of the tubes 
were then wrapped with wire having diameters ranging from 0 . 1 
to 1.6 mm to test enhancement effects. To attach the wire to 
the tubes, two different methods were used. The larger 
diameter (stiff er) wires (;» 0.5 mm) were wrapped around a rod 
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of a smaller diameter than the tube, this formed the wire into 
a tight coil. The wire was then fed onto the tube by twisting 
the coil and advancing it along the groove of the tube. The 
wire was then spot welded in place at each end of the active 
condensing surface. For the wires < 0.5 mm, one end was spot 
welded to the tube at the edge of the active surface. The 
wire was then tightly wrapped around the tube, by hand, in the 
groove and spot welded to the other end of the active surface. 
Both methods resulted in a constant wire pitch equal to the 
corrugation pitch of the tube. The wire was not physically 
attached to the active tube surface and did not act as a fin. 
The additional tubes included in Table I were manufactured for 
research conducted by O'Keefe [Ref. 18] and Mitrou [Ref. 14] . 
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Table I 



LISTING OF TUBES TESTED 
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‘Fractional coverage in this thesis is the ratio of the product of the wire diameter and 
the wire length on the active surface to the total active surface area. 



IV. EXPERIMENTAL PROCEDURES AND DATA ANALYSIS 



A. TUBE PREPARATION 

To ensure consistent results in data collection, a 
continuous film of condensate must be present over the entire 
tube surface. A chemical treatment procedure used by 
Guttendorf [Ref. 16] and several other researchers at NPS was 
used to produce filmwise condensation on the tube surface. 
Prior to installation into the test section, the tube was 
chemically treated as follows: 

1. The inside and outside of the tube were thoroughly 
cleaned with mild soap using a soft bristle brush. The 
soap was rinsed off with distilled water, then the tube 
was rinsed with acetone. Distilled water was again 
poured over the tube, to ensure the water formed a 
continuous film on the tube surface. If breaks in the 
film were observed (an indication of impurities on the 
tube surface) the cleaning process was repeated. Once a 
continuous film was present, the active surface was not 
handled, thus minimizing the possibility of 
contamination . 

2. The tube was then placed in a steam bath. 

3. A solution of 50% by weight sodium hydroxide and ethyl 
alcohol was mixed for treating the tube surface. The 
solution was kept warm to ensure a watery consistency. 

4. The solution was then applied over the tube surface, with 
a small brush, in 10 minute intervals for one hour. 
While applying the solution, the tube was rotated so the 
entire surface could be treated. If the tube was treated 
previously, the solution was only applied every 5 minutes 
over a 20 minute period. 

5. After the tube surface was treated, it was removed from 
the steam bath and rinsed with distilled water to remove 
any excess alcohol/sodium hydroxide solution. To verify 
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a continuous film will be produced in the presence of 
steam, the tube was returned to the steam bath, to ensure 
the entire surface wets easily. 

6. Without touching the active surface, the tube was removed 
from the steam bath and immediately mounted in the test 
section of the boiler apparatus. If tests were to be 
conducted with an insert, the insert was installed prior 
to securing the coolant flow piping to the test section. 



The oxide layer formed by the chemical treatment promotes 
a continuous film on the tube surface. Due to the negligible 
thickness of the oxide film, its thermal resistance was 
considered insignificant compared to the overall thermal 
resistance of the tube. 



B. SYSTEM START-UP AND SHUTDOWN PROCEDURES 

After the tube has been installed in the test section, the 
system can be brought on-line using the following procedures: 

1. The boiler section must be filled with distilled water to 
a level approximately 4 to 6 inches above the top of the 
heater elements. The boiler can be filled by attaching 
a hose from the distilled water tank to the piping at the 
boiler drain/fill valve. Opening the vent valve on the 
auxiliary condenser allows for gravity filling of the 
boiler. To drain the boiler, disconnect the distilled 
water tank hose and attach the hose leading to the drain 
in bilge area below the boiler. Then open the drain/fill 
valve and allow the water to gravity drain from the 
boiler . 

2. When the appropriate water level has been achieved, shut 
both the fill/drain valve and the vent valve on the 
auxiliary condenser. 

3. Energize the data acquisition system, computer, printer, 
and quartz thermometer power supply. Load the software 
program DRPML and check for proper operation. Then 
verify that the thermocouple outputs all correspond to 
the same ambient temperature. 
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4. Open the fill valve to the cooling water sump and adjust 
the flow rate to ensure the drain box does not overflow. 
The fill valve for the sump is located to the left of the 
boiler switchboard. 

5. Perform a leak test on the cooling water systems. 

a. Energize the cooling water supply pumps and adjust 
the rotameter flow rate from 20% to 80% while 
checking for leaks in the test section. Then secure 
the control valve and de-energize the cooling water 
supply pumps . 

b. Open the control valve for the auxiliary condenser 
and adjust the flow rate to at least 30% and check 
for leaks in the system. After completing the leak 
check, reset the flow rate to the auxiliary 
condenser to the desired position, minimum of 10%. 

6. Place the system under a vacuum. To do this, shut the 
drain valve in the plexiglass cold trap and energize the 
vacuum pump. Check the pressure gage on the vacuum pump 
to ensure proper operation and open the suction valve 
located on the side of the auxiliary condenser. After 
the vacuum pump has reduced the system pressure to below 
3 psia, shut the suction valve and secure the vacuum 
pump. 

7. With the system under a vacuum, the heaters may be 
energized . Three switches must be placed in the "ON" 
position to energize the heaters. 

a. Switch #3 on panel p5 , in the main hallway adjacent 
to room H-106. 

b. The heater load bank circuit breaker on the left 
side of the boiler heater switchboard. 

c. The condensing rig boiler power switch, on the front 
of the boiler heater switchboard. 

When the heaters are energized, the power level should be 
adjusted to 50 Volts, using the silicon control rectifier 
on the boiler control switchboard. If the system is below 
2 psia, the power level should initially set to 40 Volts 
to minimize vibrational shock to the system caused from 
vapor bubble formation and collapse. As the system warms 
up, the power can be increased, in increments of 10 
Volts, until the desired setting is reached. 
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8. After the system pressure rises above 4 psia, non- 
condensible gases need to be purged from the system. 
This can be accomplished by following the procedures 
outlined in step 6. Throughout the purging procedure, 
the horizontal tube coolant flow should remain secured, 
while the flow rate to the auxiliary condenser can be 
adjusted until all the gases have been evacuated. When 
the auxiliary condenser is warm to the touch everywhere, 
and condensate has formed in the plexiglass container, 
most of the non- condensible gases have been removed. The 
initial purging of the system may take between 15 and 30 
minutes. This process should be repeated every few hours 
during extended operation. 

9. To establish filmwise condensation on the tube being 
tested, the following steps should be performed: 

a. Allow the system vapor temperature to increase to 
=90 °C (at least 3800 microvolts on channel 40) . 

b. Increase the flow rate in the auxiliary condenser 
to 50% or 60% and allow the vapor temperature to 
reduce to =75 °C (approximately 3200 microvolts) . 

c. Secure coolant flow to the auxiliary condenser and 
allow the vapor temperature to increase to about 85 
°C (=3700 microvolts) . This forms a steam blanket 
over the tube. 

d. Initiate coolant flow of 80% in the single 

horizontal tube. 

e. Adjust the auxiliary condenser flow rate to maintain 
the desired temperature and pressure for data 
collection . 

If, after performing steps 9 a. through 9 e., some 
dropwise condensation persists, repeat step 9 again. If 
dropwise condensation continues, the tube should be 
removed, cleaned and retreated in accordance with the 
procedures listed in section IV. A. A cutout section in 
the insulation above the test section allows for a light 
to be shown on the upper surface of the tube in order to 
check for a continuous film on the tube surface. 

10. Run the software program DRPML by pressing the "RUN" 
key on the keyboard. The program will prompt you 
with questions for the necessary information 
required to run as listed by O'Keefe [Ref. 18] and 
in Appendix C. 
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11. Ensure the system has been operating at steady- state 
conditions for at least 30 minutes prior to 
continuing past the question "Enter the flowmeter 
reading" . 

12. Vacuum runs are conducted at a heater setting of 90 
Volts and 1980 ± 10 microvolts on channel 40. This 
corresponds to T Sat c=48 C » and a vapor velocity of ■== 
2 m/ s . 

13 . The atmospheric runs are conducted at a heater 

setting of 175 Volts and 4280 ± 10 microvolts on 

channel 40. This corresponds to T Sat =100°C, and a 
vapor velocity of <=1 m/s. Special care must be used 
when operating at atmospheric pressure to ensure 
over- pressurization and rupture does not occur . 

14. When taking readings, double check the flowmeter 
reading prior to accepting any data point. If both 
vacuum and atmospheric pressure runs are to be 
conducted in the same day, the vacuum data should be 
collected first. This eliminates the long cool down 
time required after an atmospheric pressure run. 

After data has been collected, the system should be 
secured utilizing the following procedures: 

1. Secure power to the heating elements. Turn off the 
switches in the reverse order listed in step 7 of the 
start-up procedures. 

2. Secure coolant flow to the auxiliary condenser. If the 
system is to remain at vacuum pressure until the next 
data run, then the auxiliary condenser can be used in 
assisting to cool the system down, provided the same tube 
will be used for the following test. 

3. Secure the coolant flow through the single horizontal 
tube by securing the coolant pumps . 

4. Secure the water flow to the coolant sump tank. 

5. To return the system to atmospheric pressure, slowly open 
the vent valve on the auxiliary condenser. To avoid 
contaminating the system, keep any foreign material away 
from the vicinity of the vent valve during this 
procedure . 

6. If an emergency should arise, such as an 
over-pressurization or breakage, ensure the heater power 
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is secured first. Then, let the system cool down prior 
to checking for damage. 

C. DATA COLLECTION PROCEDURE 

After the tube has been installed and the system has been 
operating for a minimum of 30 minutes, data collection may 
begin. Each run consisted of data points taken at coolant 
flow rate settings from 80% to 20% and back up to 80%, in 10% 
increments. This allowed for each data point to be checked 
for repeatability during the run. Before a data point was 
accepted, sample readings were taken to ensure the coolant 
temperature difference, the overall heat transferred, and the 
saturation temperature were stable. The temperature 

difference between the inlet and outlet coolant flow must be 
stabilized (+0.01 Q C) prior to recording data . After setting 
the coolant flow rate, the aforementioned parameters required 
between 5 and 15 minutes to stabilize. 

D. DATA REDUCTION PROCEDURE 

For a slice of the condenser tube, the overall thermal 
resistance from vapor to coolant can be represented by the sum 
of the coolant side resistance (R±) , the fouling resistance 
(R f ) , the wall resistance (R w ) , and the vapor side resistance 
(R Q ) , as shown in equation (4.1). 

R Cotal = R i +R f +R w +R o ( 4 ‘ 1 > 
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Only clean tubes were tested so the fouling resistance was 
set to zero, therefore: 

R tocal = R i +R w+ R o < 4 * 2 > 

Being convective in nature, the vapor and coolant side 
resistances can be expressed as follows: 



R 



o 



1 

hoA 0 



(4.3) 



R i = 



h i A 1 



(4.4) 



The wall resistance term comes from the general radial 
conduction equation for a concentric cylinder: 



= ^DjDj) 

2k Lk m 



(4.5) 



The total thermal resistance can be expressed in terms of the 
overall heat transfer coefficient (U Q ) and the effective 
outside condensing area (A Q ) , as follows: 



R 



total 



l 

u<Ao 



(4.6) 



The effective outside condensing area of the tube can be 
calculated using the following equation: 



A 0 = 7 iD 0 L (4.7) 

For the inside effective area, two different conditions 
are present. The area corresponding to the outside effective 
area, where heat transfer takes place radially, and the inlet 
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and outlet portions of the tube, where heat removal occurs 
axially. To account for this 'fin' effect, the extended fin 
assumption and associated fin efficiencies were used in the 
following expression: 

Aj = nD i (L+L 1 ri 1 +L 2 ri 2 ) (4.8) 

Substituting equations (4.3), (4.4), and (4.6) into 

equation (4.2) yields the following: 



1 

Wo 



h 1 A i 



+ R, 



l 

hoAo 



(4.9) 



From the data collected, the total heat transfer rate 
across the tube can be calculated directly from an energy 
balance using the coolant mass flow rate and the coolant 
temperature difference (from inlet to outlet) through the 
tube : 



Q = mc p {T 2 -T 1 ) (4.10) 

The total heat transfer rate can also be expressed in terms of 
the overall heat transfer coefficient, the effective outside 
condensing area, and the log mean temperature difference 
(LMTD) : 



Q = U 0 A 0 ( LMTD) 



where : 



LMTD 



(r 2 -Ti) 



In 



T -T 

_ sa c ■ L i 

T -T 
1 sac 1 2 



(4.11) 



(4.12) 
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The measurements for LMTD were obtained in the following 



manner : 



• T 1 was measured from a quartz thermometer installed at the 
inlet section of the coolant flow for the condenser tube. 

• T 2 was measured with a quartz thermometer inserted 
downstream of the mixing chamber at the outlet section of 
the coolant flow for the condenser tube. A correction 
factor was applied to T 2 in order to compensate for 
viscous heating of the coolant through the tube. The 
corrections, which depend on the inside geometry of the 
tube, the coolant flow velocity, and the type of insert, 
are listed in Appendix A. 

• T sat was measured using the vapor thermocouple (channel 
40) . 



Substituting equation (4.10) into equation (4.11) provides 
determination of the overall heat transfer coefficient 
directly from experimental measurements: 



U„ 



ihc p (T 2 -T 1 ) 
AJLMTD) 



( 4 . 13 ) 



With A Q , A it R w , and U Q all known quantities, the only 
remaining unknowns in equation (4.9) are h Q and h i? the 
outside and inside heat transfer coefficients. The modified 
Wilson plot technique was used to calculate these 
coef f icients . 



E. MODIFIED WILSON PLOT TECHNIQUE 

In 1915, Wilson [Ref. 63], developed a method for 
indirectly determining the inside and outside thermal 
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resistance from an overall resistance. 



Wilson's method 



requires the heat flux of the system to remain constant . 
While this application proves valid in theory, difficulties 
arise when trying to maintain a constant heat flux while 
varying the cooling water velocity, Marto [Ref. 64] . 

Khartabil et al . [Ref. 65] provide a comprehensive listing 
of modifications to the original technique developed by 
Wilson. For a detailed development of the modified Wilson 
plot technique, the reader is referred to Marto [Ref. 64] and 
Marto and Nunn [Ref. 66] . 

Using the modified Wilson plot method to determine the 
inside and outside heat transfer coefficients requires one to 
choose a general form for each correlation. The program DRPML 
enables the user to select from several options. The inside 
heat transfer coefficient correlation can be written in 
simplified form as follows: 

^ = CVQ (4.14) 

where Q varies with the particular correlation used. Using 
the Sieder-Tate correlation [Ref. 49] (equation (2.19)): 

Q = — Re x Pr 1/3 ( |i r /|i J °- 14 (4.15) 

Di 

where the program allows the Reynolds number exponent (x) to 
be varied. Using the Sleicher-Rouse correlation [Ref. 50] 
(equation (2.21) ) : 
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(4.16) 



Q 



K 

1 



5 + 0 . 015 Re f c Pr w d ) 



where c and d are defined in equation (2.22) and (2.23) 
respectively. Using the Petukhov- Popov correlation [Ref. 51] 
(equation (2.24)): 



Q 




( e / 8 ) Re Pr_ 

+ K 2 (e/8) 1/2 (Pr 2/3 -l) 



(4.17) 



where e, , and K 2 are defined in equation (2.25), (2.26), 

and (2.27) respectively. 

Similarly, the outside heat transfer coefficient 
correlation may be written in simplified form: 

h a = a Z (4.18) 

where Z can be selected from equation (2.1) (Nusselt's theory 
[Ref . 23] ) : 



Z = 



*1 9 P f 2 V 

(T sat~T wo ) 



(4.19) 



or from equation (2.62) (Fujii's correlation [Ref. 29]): 

Z = ^ F 1/5 Re 3/2 (4.20) 

D 0 241 

F and are defined in equation (2.4) and (2.5) 

respectively . 
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Substituting equation (4.14) and (4.18) into equation 



(4.9) yields: 



1 _ 1 1 
~U^A 0 ~ C i QA i + Rw + 

which can be written as: 



( 4 . 21 ) 



1 


1 


A.Z] 


~o A 


• z = — • 

Ci 


QA X 



JL 

a 



By letting: 



Y = 




• Z 



( 4 . 22 ) 



( 4 . 23 ) 



X = 



AgZ 

A t Q 



( 4 . 24 ) 






Ci 



and : 



b = 



1 

a 



( 4 . 25 ) 



( 4 . 26 ) 



Equation (4.22) can be written in simple linear form as: 

Y = m x X + b ( 4 . 27 ) 

Due to the temperature dependence of the parameters Q and 
Z, the program utilizes an iterative procedure in the solution 
for each data point. A least - squares fit of equation (4.27) 
then determines and a. To calculate h Q , the known value of 
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a can be used directly in equation (4.18) ; alternatively, with 
the known value of C^, h i can be calculated from equation 
(4.14), and applied to equation (4.9) in the following form: 




-i 



h 



( 4 . 28 ) 



O 



F. ENHANCEMENT RATIO 

The enhancement ratios in this thesis are used to express 
a comparison of enhanced tube results to smooth tube results, 
or to compare wire-wrapped corrugated tube results to plain 
corrugated tube results. 

From Nusselt's theory, heat flux can be written as: 



q = a (T sac -T wo ) 3/4 



( 4 . 29 ) 



where : 




1/4 



( 4 . 30 ) 



Heat flux can also be represented by: 



or = K (T sac -T wo ) 



( 4 . 31 ) 



Equating equations (4.29) and (4.31) yields: 
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h = a (T -T ) ~ 1/4 

11 o a ' 1 sac 1 wo' 



( 4 . 32 ) 



Taking the ratio of the outside heat transfer coefficient 
from an enhanced tube, (h oe ) , to a smooth tube (h os ) , for the 
same (T sac -T wo ) yields: 



h 



e T 



oe _ 



h. 



a. 

a. 



a. 
a „ 



( 4 . 33 ) 



Similarly, the enhancement ratio for a wire -wrapped corrugated 
tube to a plain corrugated tube at the same (T sat -T wo ) can be 
represented by: 



-Tc 



^owc _ ^ wc 



CL, 



( 4 . 34 ) 



where the subscripts wc and c refer to wire-wrapped corrugated 
and plain corrugated tubes respectively. 

For the purpose of this study all enhancement results for 
wire-wrapped corrugated tubes are based on e Tc and wire- 
wrapped smooth tube enhancements are based on e T . 
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V. RESULTS AND DISCUSSION 



A. HEAT TRANSFER CORRELATIONS 

Lorenz et al . [Ref. 52] conducted a thorough investigation 
of several inside heat transfer correlations to determine 
which correlation most accurately predicted their data for 
turbulent flow of water in smooth tubes. Their findings 
showed that the Petukov- Popov and Sleicher- Rouse correlations 
(equations (2.24) and (2.21) respectively) were the most 
accurate in predicting the inside heat transfer coefficient 
over the range of Prandtl numbers from 6.0 to 11.6. 
Consequently, O'Keefe [Ref. 18] made a comparison of his 
results obtained using both correlations to process data, 
which was collected from the same experimental apparatus used 
in this study. His results showed good agreement between the 
outside coefficients obtained for smooth copper and smooth 
titanium tubes, when the two correlations were used with the 
modified Wilson plot data reduction technique. Based on his 
findings, O'Keefe [Ref. 18] processed most of his data using 
the Petukhov- Popov correlation. All data collected in this 
study were processed using the Petukov- Popov correlation so a 
direct comparison could be made with the results obtained by 
O'Keefe fRef. 181 . 
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The data collection and reduction program, DRPML, listed 
in Appendix C, allows for the selection from two possible 
outside heat transfer correlations, the Nusselt theory or the 
Fujii correlation (equations (2.1) and (2.6) respectively). 
Figure 5 shows the reduced results for h 0 obtained from using 
these two equations for a vacuum pressure run on a plain MHT 
Korodense tube. It can be seen that there is good agreement 
between the two equations. This indicates that the leading 
coefficient, oi, obtained from the modified Wilson plot method 
using the Nusselt theory compensates for any vapor velocity 
effect inherent in the Fujii correlation. O'Keefe [Ref. 18] 
processed all of his data using the Nusselt correlation; so 
for comparison purposes, all further calculations of the 
outside heat transfer coefficient in the present work have 
also been made using the Nusselt theory. 

B. VERIFICATION AND REPEATABILITY OF RESULTS 

Several runs were made using the same tubes tested by 
previous researchers to verify system operation and data 
collection procedures. The five tubes that were tested 
included two smooth titanium tubes and two LPD Korodense 
titanium tubes tested by O'Keefe [Ref. 18]; of these, one of 
each was wrapped with a 0.5 mm diameter titanium wire. The 
fifth tube was a smooth copper tube wrapped with a 0.5 mm 
diameter titanium wire from Mitrou's work [Ref. 14] . Figures 
6 and 7 are comparisons of the outside heat transfer 



49 



Ho vs. Tcf 

MHT Korodense Tube Under Vacuum 




CD 



CD 



CD 



(N-Z.w/M*) oh 



Figure 5. h Q vs . AT cf for a Plain MHT Korodense Tube Under 
Vacuum 
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T cl (K) 



Ho vs. Tcf 

Smooth Titanium Tube Under Vacuum 




(>I-2.uj/m>I) °H 



Figure 6 . 



h Q vs. AT cf for a Smooch Titanium Tube Under Vacuum 
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Tcf (K) 



Ho vs. Tcf 

Smooth Titanium Tube Atmospheric Pressure 




(N-3.UJ/AM) OH 



Figure 7. h Q vs . AT cf for a Smooth Titanium Tube at 
Atmospheric Pressure 



52 



Tcf (K) 



coefficient for a vacuum and an atmospheric pressure run 
respectively for data collected in this study and the study 
conducted by O'Keefe [Ref. 18] for the smooth titanium tube. 
The figures show good agreement for the two sets of data. The 
separation in the data sets seen in Figure 6 is due to a 
higher coolant water inlet temperature when O'Keefe conducted 
his runs. Tables II through VII contain a listing of all the 
data collected by O'Keefe, Mitrou and from the present study 
for comparison purposes in this study. The researcher 
initials are as follows: (0)=0'Keefe [Ref. 18], (L)=Long, and 
(M)=Mitrou [Ref. 14]. The average variation for a, the 
leading coefficient for the outside heat transfer equation, 
was less than 3.3%, with 90% of the results differing by less 
than 2%. The C i values obtained in this study differed by as 
much as 8.5% compared with values obtained by previous 
researchers. This larger difference is possibly due to the 
uncertainty as to which HEATEX insert was used in O'Keefe's 
experiments. Four different inserts were available and each 
had a slightly different length. 

C. ANALYSIS OF CORRUGATED TUBE RESULTS 

To establish a base line for the enhancement ratios of the 
wire-wrapped MHT Korodense tubes, an initial series of runs 
were conducted to determine the heat transfer characteristics 
of the plain MHT Korodense tube. The data obtained are 
presented in Table VIII. Figure 8 shows a near unity 
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Table II. SMOOTH TITANIUM TUBE AT ATMOSPHERIC PRESSURE 







HEATEX 


Insert 


Data Run/Researcher 


C i 


Of 


F0NMAHT1 


(0) 


2 .31 


0 . 71 


F0NMAHT2 


(0) 


2 .38 


0.77 1 


F0NMAHT3 


(0) 


2.39 


0 . 76 


F0NMAHT4 


(0) 


2 .41 


0.75 


F0NMAHT5 


(0) 


2.20 


0 . 77 


F0NMAHT6 


(0) 


2 . 54 


0.76 


F0NMAHT7 


(0) 


2.51 


0.75 


LSNTHA1 


(L) 


2.90 


0.75 






No 


Insert 


F0NMANT1 


(0) 


1.21 


0.75 


F0NMANT2 


(0) 


1.18 


0.76 


F0NMANT3 


(0) 


1 . 17 


0 . 77 


F0NMANT4 


(0) 


1 . 18 


0 . 76 


F0NMANT5 


(0) 


1.24 


0 . 79 


LSNTNA1 


(L) 


1.17 


0.82 



enhancement of the outside heat transfer coefficient for the 

MHT Korodense tube (having a groove pitch of 8 mm) over smooth 

titanium tube results. This is consistent with the findings 

of Rabas [Ref. 41] where, in a comparison of an extensive 

collection of data on spirally indented tubes he found: 

. . . that there is no appreciable enhancement for 
condensation on roped tubes for groove pitches greater 
than 4 mm. 
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Table III. SMOOTH TITANIUM TUBE AT VACUUM PRESSURE 







HEATEX 


Insert 


Data Run/Researcher 


C, 


Ct 


F0NMVHT3 


(0) 


2 . 55 


0.15 


F0NMVHT4 


(0) 


2.29 


0 . 79 


F0NMVHT5 


(0) 


2.42 


0 . 19 


LSNTHV1 


(L) 


2 . 87 


0.78 


LSNTHV2 


(L) 


2 . 11 


0 . 19 


LSNTHV3 


(L) 


2.43 


0 . 79 


LSNTHV4 


(L) 


2 . 50 


0.7* 


No Insert 


F0NMVNT2 


(0) 


1.03 


0.82 


F0NMVNT3 


(0) 


1.08 


0 . 79 


F0NMVNT4 


(0) 


1.07 


0.85 


F0NMVNT5 


(0) 


1 . 11 


0.82 


LSNTNV1 


(L) 


1 . 13 


0.83 


LSNTNV2 


(L) 


1.13 


0.84 


|| LSNTNV4 


(L) 


1.11 


0 . 83 



A comparison of the Ci values for the no insert runs on 
the smooth tube (Tables II and III) and LPD tube (Table V) 
indicate increases in h^ from 73 to 74% for vacuum and 
atmospheric pressure runs respectively. As discussed earlier, 
Wolverine [Ref. 2] reported a 25% increase in the inside heat 
transfer coefficient for the MHT Korodense tube when compared 
to the LPD Korodense tube. Comparison of the C i values for the 
no insert runs on the MHT tube (Table VIII) and the LPD tube 
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Table IV. 0 . 5 mm WIRE -WRAPPED SMOOTH TITANIUM TUBE WITH 
HEATEX INSERT 



Atmospheric Pressure 


Data Run/Researcher 


C i 


a 


FONMAH6T2 (0) 


2 .39 


0.99 


FONMAH5T3 (0) 


2 .37 


1 . 01 


LSW1THA1 (L) 


2 . 67 


1 . 04 


Vacuum Pressure 


F0NMVH6T1 (0) 


2 . 16 


0 . 95 


FONMVH6T2 (0) 


2.21 


0.95 


LSW1THV1 (L) 


2.59 


0.94 



(Table V) indicate increases in h i comparable to those 
reported by Wolverine [Ref. 2], 21% for the atmospheric 
pressure runs and 25% for the vacuum runs. This indicates 
that as the groove depth increases from the smooth tube (zero) 
to the LPD tube value (<=0.33 mm) most of the increase in the 
inside heat transfer coefficient due to corrugation has 
already been realized. The additional groove depth of the MHT 
tube provides only a 21 to 25% further increase of the inside 
heat transfer above the LPD results (MHT tube groove depth 
<=0.5mm), while the pressure drop doubles [Ref. 2]. 

Figures 9 and 10 show the effect of using the HEATEX 
insert on the overall heat transfer coefficient (U 0 ) • The 
general trend for U 0 to increase with increasing coolant flow 
velocity is due to improved coolant mixing and a lower mean 
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Table V. PLAIN LPD KORODENSE TUBE 



Atmospheric Pressure HEATEX 


Insert 


Data Run/Researcher 




Of 


FONMAHLT2 


(0) 


2 . 90 


0.90 


FONMAHLT3 


(0) 


2 . 89 


0 . 92 


LLNTHA1 


(L) 


3 .11 


0.95 




Atmospheric Pressure 


No Insert 


FONMANLT2 


(0) 


2.06 


0.92 


FONMANLT3 


(0) 


1.99 


0.94 


LLNTNA2 


(L) 


2 . 17 


0.95 


Vacuum Pressure HEATEX Insert 


FONMVHLT1 


(0) 


2 . 72 


0.94 


F0NMVHLT2 


(0) 


2 .67 


0.95 


LLNTHV2 


(L) 


3 . 02 


0.97 




Vacuum Pressure No 


Insert 


FONMVNLT2 


(0) 


1.87 


0.95 


FONMVNLT3 


(0) 


1.86 


0.96 


LLNTNV3 


(L) 


1.99 


0.96 



bulk fluid temperature. Increases range from 3% at a coolant 
flow rate of 3.83 m/s up to 9% at a flow rate of 1.02 m/s. 
This increase is attributed to the greater turbulence of the 
coolant flow due to the presence of the insert and can be seen 
in the higher values for the HEATEX runs for every tube 
tested. A comparison of O'Keefe's data [Ref. 18] for smooth 
titanium and LPD Korodense titanium tubes, tested with and 
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Table VI. 0 . 5 mm WIRE - WRAPPED LPD KORODENSE TUBE WITH 
HEATEX INSERT 



Atmospheric Pressure 


Data Run/Researcher 




a 


FONMAHL3T3 (0) 


2 . 84 


0.97 


LLW05THA1 (L) 


3.03 


0.98 


LLW05THA2 (L) 


3.06 


o 

CD 


Vacuum Pressure 


F0NMVHL3T1 (0) 


2.61 


0.99 


LLW05THV1 (L) 


2 . 78 


0.98 


LLW05THV2 (L) 


2 . 78 


0.98 



without a HEATEX insert indicates a 20% increase in U 0 for the 
smooth tube and a 7.9% increase in U 0 for the LPD tube when 
tested at a coolant velocity of »=2 . 5 m/s. For a coolant flow 
rate of =2 . 5 m/s the present study shows a 7.2% increase in U 0 
for the MHT tube when a HEATEX insert is used. This trend 
indicates that a more severe corrugation groove depth (zero 
for the smooth tube, 0.33 mm for the LPD tube and 0.5 mm for 
the MHT tube) results in a smaller increase in U 0 with the 
addition of an insert. The grooves present in the MHT and LPD 
tubes cause greater turbulence in the coolant flow than the 
inside walls of the smooth tube. Therefore, the addition of 
an insert for a corrugated tube (which already has turbulent 
flow) has a less pronounced effect on U 0 compared to the 
effect realized for a smooth tube. 
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Table VII. 0 . 5 mm WIRE-WRAPPED SMOOTH COPPER TUBE WITH 
HEATEX INSERT 







Atmospheric 


Pressure 


Data Run/Researcher 


C i 


Of 


FONMAH71C1 


(0) 


3 . 07 


1.50 


LSW71CHA1 


(L) 


3 . 15 


1.28 


S71A305* 


(M) 


2 . 69 


1.50 


S71A314* 


(M) 


2 . 69 


1.50 


Vacuum Pressure 


F0NMVH71C2 


(0) 


2 . 76 


1.28 


LSW71CHV1 


(L) 


2.94 


1.23 


S71V296* 


(M) 


2.50 


1.27 



Runs conducted with wire -wrapped insert 



D. ANALYSIS OF WIRE-WRAPPED CORRUGATED TUBE RESULTS 

Seven different titanium wire diameters ranging from 0.1 
mm to 1.6 mm were used to wrap separate MHT Korodense tubes. 
Tests were conducted at atmospheric and vacuum conditions and 
the results were compared to the plain MHT Korodense tube data 
to determine the enhancement effects attributable to wire- 
wrapping. Wire -wrapped smooth tube data obtained by previous 
researchers were also used for comparison to establish if 
similar trends exist. 

The enhancement due to wire-wrapping can be seen directly 
in Figure 11 where the averaged U 0 (averaged for all data runs 
at each flow rate) are plotted for all the wire diameters. 
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Table VIII. PLAIN MHT KORODENSE, TUBE 



Atmospheric Pressure HEATEX Insert 


Data Run 


c,- 


(2 


LMNTHA7 


3 . 34 


0.79 


LMNTHA8 


3.31 


0.79 


LMNTHA9 


3.32 


0.80 


Atmospheric Pressure No Insert 


LMNTNA7 


2.49 


0.80 


LMNTNA8 


2.47 


0.80 


LMNTNA9 


2 .56 


0.80 


Vacuum Pressure HEATEX Insert 


LMNTHV7 


3 . 07 


0.78 


LMNTHV8 


3 . 01 


0.79 


LMNTHV9 


3.06 


0.78 


Vacuum Pressure No Insert 


LMNTNV7 


2.36 


0.77 


LMNTNV8 


2 . 39 


0.76 


LMNTNV9 


2.43 


0 . 77 



Here the 0.4 mm wire-wrap appears to give the best heat 
transfer with an enhancement of about 3%. Similar trends for 
the outside heat transfer coefficient can also be seen in 
Figures 12 and 13 . The plain MHT Korodense tube results are 
included to show the relative enhancements for each size of 
wire. Figure 12 shows an increasing trend with decrease in 
wire diameter starting from the 1.6 mm wire -wrap (P w /D w =5) and 
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MHT Korodense Tube and Smooth Titanium Tube Under Vacuum 




(>i-2.uj/M>|) oh 



Figure 8. h Q vs . AT cf for a Smooth Titanium and a Plain MHT 
Korodense Tubes Under Vacuum 
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Tct (K) 



Uo vs. V w 

Plain MHT Korodense Tube Atmospheric Pressure 
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Figure 9. U 0 vs . V w for a Plain MHT Korodense Tube at 
Atmospheric Pressure 
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Coolant Flow Velocity (m/s) 



Uo vs. Vw 

Plain MHT Korodense Tube Under Vacuum 
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Figure 10. U Q vs. V w for a Plain MHT Korodense Tube Under 
Vacuum 
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Coolant Flow Velocity (m/s) 



Uo vs. Vw 

Wire-Wrapped MHT Korodense Tube Under Vacuum with HEATEX 

Insert 
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Figure 11. U Q vs. V w for Wire-Wrapped MHT Korodense Tubes 
Under Vacuum with HEATEX Insert 
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Flow Velocity (m/s) 



Ho vs. Tcf 
Under Vacuum 
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Figure 12. h Q vs. T cf for Wire-Wrapped MHT Korodense Tubes 
Under Vacuum 
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Tcf (K) 



an enhancement ratio of 0.88 to the 0.4 mm wire -wrap 
(P w /D w =20) with an enhancement ratio of 1.08. Figure 13 
completes the data and shows the trend returning to the plain 
MHT tube results, starting with the 0.4 mm wire- wrap and 
ending with the 0.1 mm wire-wrap (P w /D w =80). Two separate 
figures were used to help eliminate the clutter from the 
overlapping data. Table IX contains the averaged C i and a 
values for all the conditions tested (vacuum and atmospheric 
pressure, with and without the HEATEX insert) for the plain 
MHT tube and all the wire -wrapped MHT tubes. Two to four 
tests were conducted at each condition to ensure repeatability 
of the data. 

The wire-wrapped MHT tube data were also checked to see if 
there was a relationship between the fractional coverage and 
the pitch to wire diameter ratio to the enhancement over the 
plain MHT tube (E Tc ) . Figure 14 shows the relationship 
between the enhancement ratio and the pitch to wire diameter 
ratio and Figure 15 shows enhancement ratio versus fractional 
coverage. Data obtained by O'Keefe [Ref. 18] for wire-wrapped 
LPD Korodense tubes are also included in the figures. 
However, the number of his tests are limited to only 3 wire 
diameters (0.5mm, 1.0mm and 1.6mm) . The data clearly indicate 
that wrapping with a large diameter wire causes a degradation 
in performance. As the wire diameter is decreased, the 
enhancement passes through a maximum of about 8%, 
corresponding to a P w /D w «=20 and a fractional coverage of =5%. 
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Table IX. WIRE-WRAPPED AND PLAIN MHT KORODENSE TUBES 



Atmospheric Pressure 


Tube 


No Insert 


Heatex Insert 


c i 


Of 


Ch 


Of 


Plain 


2.51 


o 

CD 


3 . 32 


0 . 79 


1.6mm wire -wrap 


2 .31 


0 . 72 


2.92 


0.73 


1.0mm wire-wrap 


2 .35 


0 . 78 


3.17 


0 . 79 


0.5mm wire -wrap 


2 . 50 


0 .85 


3 . 19 


0 . 85 


II . 0mm wire-wrap 


2 . 53 


0 . 84 


3 .38 


0 . 85 


0.3mm wire -wrap 


2.52 


0 . 84 


3.27 


0 . 85 


a . imm wi re - wrap 


2 . 53 


0 . 82 


3 .31 


0 . 82 


0.4mm wire-wrap 


2 . 50 


0.81 


3 .36 


0.82 


Vacuum Pressure 


Plain 


2.99 


0 . 77 


3.05 


0 . 79 


1.6mm wire -wrap 


2 . 17 


0 . 67 


2 . 84 


0.69 


1 . 0mm wi re - wrap 


2.20 


0.12 


2.99 


0 . 73 


0.5mm wire -wrap 


2.14 


0.82 


3.05 


0 .81 


0.4mm wxre-wrap 


2.31 


0.83 


3 . 03 


0.84 


0.3mm wire -wrap 


2.31 


0.82 


3 . 14 


0 . *3 


0 . 2mm wire - wrap 


2 . 39 


0 . 81 


3 . 16 


0 . 80 


0.1mm wire -wrap 


2.42 


0 . 78 


3.23 


0 . 80 



As the wire diameter is further decreased, the enhancement 
ratio returns to unity, i.e. back to the plain tube 
performance. The results from O'Keefe [Ref. 18] for the LPD 
Korodense tube definitely show the same trends as the MHT 
tube, but the results are somewhat inconclusive due to lack of 
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Enhancement Ratio vs. Pitch to Diameter Ratio 
Wire-Wrapped MHT and LPD Korodense Tubes 
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Figure 14. Enhancement Ratio vs. P w /D w for Wire-Wrapped 
Korodense Tubes 
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Pw/Dw 



Enhancement vs. Fractional Coverage 
Wire-Wrapped MHT and LPD Korodense Tubes 



co 

o 






E 

cd 



2 B 

> cd 



o 

cd 

> 



a 

a. 



a 

CL 




i 




! 



a> 

o 



o i j b y )U9tu9ouBi|ug 



CM 

O 



O 

O 



GO 

O 



Figure 15. Enhancement Ratio vs. Fractional Coverage for 
Wire -Wrapped Korodense Tubes 
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Fractional Coverage 



data. His smallest wire diameter was only 0.5 mm. 

The enhancement results obtained for the wire-wrapped MHT 
Korodense tube were also compared to a wire -wrapped smooth 
tube. Marto et al . [Ref. 31] reported an optimum wire pitch 
to diameter ratio ranging from 5 to 7 for wire-wrapped copper 
tubes. Sethumadhavan and Rao [Ref. 34] realized a maximum 
enhancement at a fractional coverage of 21%, corresponding to 
a P w /D w =5 . This value for P w /D w is much smaller than the 
optimum pitch to diameter ratio of 20 obtained in this study. 

To further understand this relationship, data were 
compiled for wire -wrapped smooth tubes studied at NPS and 
compared with the current data in Figures 16 and 17. Both 
sets of data were plotted with constant pitch curves for a 
better comparison with the constant pitch results of the 
corrugated tubes. It can be seen that the wire-wrapped smooth 
copper and wire -wrapped smooth titanium tubes with P w =4 mm 
have roughly the same slope. The same trend can be seen for 
the wire-wrapped smooth titanium and the wire-wrapped LPD and 
MHT tubes with a pitch of about 7 mm. Due to the limited data 
it is difficult to predict the optimal fractional coverage of 
the wire- wrapped smooth tubes for a constant wire pitch. 
However, it does appear to be less than the 21% obtained by 
Sethumadhavan and Rao [Ref. 34] . It can also be seen that the 
magnitude of the enhancements tend to vary with tube material 
and tube geometry. Mitrou's data [Ref. 14] for wire-wrapped 
smooth copper tubes show results 45 to 55% higher than the 
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Enhancement Ratio vs. Pitch to Diameter Ratio 
Wire-Wrapped Smooth and Korodense Tubes 



o 

CD 




13 

O h- 



■ 

^ m s 

■ ■ 

h a. a. 



-c .cr si 



o 


o 


o 


o 


o 


o 


B 


B 


e 


CO 


CO 


CO 



i ! t 



Q 

CL 




O 

GO 



| 






o ij b y lueueouBijug 



Figure 16. Enhancement Ratio vs. P w /D w for Wire-Wrapped 
Korodense and Wire-Wrapped Smooth Tubes, NPS Data 
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Enhancement Ratio vs. Fractional Coverage 
Wire-Wrapped Smooth and Korodense Tubes 
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Figure 17. Enhancement Ratio vs. Fractional Coverage for 
Wire -Wrapped Korodense and Wire-Wrapped Smooth Tubes, NPS Data 
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Fractional Coverage 



wire- wrapped smooth titanium tube results for an equal pitch 
(P w =4 mm) obtained by O'Keefe [Ref. 18] . This indicates that 
tubes made of materials with high thermal conductivities (such 
as copper) may provide additional enhancement characteristics. 
The effect of tube geometry can also be seen in the comparison 
of the wire-wrapped smooth titanium tube and the wire-wrapped 
LPD and MHT tubes (all with pitches =7 mm) . The wire-wrapped 
smooth titanium tube achieved an enhancement =10% higher than 
the corrugated tubes of similar pitch. Figures 18 and 19 are 
identical to Figures 16 and 17 with the addition of the data 
obtained by Sethumadhavan and Rao [Ref. 34] . The tube 
material used in their tests is not known, but appears to have 
a thermal conductivity between copper and titanium (such as 
copper/nickel or aluminum) . The most noticeable differences 
in their data and the data from research at NPS are the 
fractional coverage of the apparent maximum enhancement and 
the opposite slope of the constant pitch curves. This is most 
likely due the wire material used for wrapping their tubes. 
They used copper wire and all the NPS tubes were manufactured 
with titanium wire . The use of a highly conductive material 
may have produced more of a 'fin' effect with the wire and 
caused the peak enhancement to shift. Further study on the 
effect of the wire material is required to determine the 
effects . 

Possible effects on the inside heat transfer coefficient 
due to wire- wrapping both smooth and corrugated tubes were 
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Enhancement Ratio vs. Pitch to Diameter Ratio 
Wire-Wrapped Smooth and Korodense Tubes 
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Enhancement Ratio vs. P w /D w for Wire- Wrapped 
and Wire -Wrapped Smooth Tubes 
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Enhancement Ratio vs. Fractional Coverage 
Wire-Wrapped Smooth and Korodense Tubes 
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Figure 19. Enhancement Ratio vs. Fractional Coverage for 
Wire -Wrapped Korodense and Wire -Wrapped Smooth Tubes 
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Fractional Coverage 



also investigated. Figure 20 shows a possible decreasing 
trend in Ct with increasing fractional coverage. However, the 
magnitude of the decrease is minimal and the overall variation 
in is within the uncertainty of the plain tube results. 

In summary, the enhancement of the outside heat transfer 
coefficient for a titanium MHT Korodense tube over a smooth 
titanium tube was close to unity. The primary enhancement in 
the overall heat transfer coefficient was found to be due to 
an increase in the inside heat transfer coefficient, 21 to 26% 
above the LPD Korodense tube coefficient and 110 to 117% above 
the smooth tube coefficient. The effect of wire- wrapping 
corrugated tubes provides minimal additional enhancement to 
the outside heat transfer coefficient, a maximum of <=8% with 
a P w /D w =20 and a corresponding fractional coverage of 5%. 
However, the addition of wire-wrap into the grooves of a 
corrugated tube may have greater benefit when used in a tube 
bundle arrangement due to improved condensate drainage. 
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Figure 20 . C ^ vs. Fractional Coverage for Wire-Wrapped 

Korodense and Wire -Wrapped Smooth Tubes at Atmospheric 
Pressure 
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Fractional Coverage 



VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. Enhancement of the outside heat transfer coefficient due 
to the corrugation of the MHT Korodense titanium tube is 
near unity when compared to smooth tube results. 

2. The corrugation of the MHT Korodense tube produces 
enhancements in the inside heat transfer coefficient 
from 110% to 117% above smooth tube results, and 21% to 
26% above LPD Korodense tube results. 

3. Wire- wrapping of MHT Korodense tubes yields minimal 
additional enhancement, «8%, when used in a single tube 
arrangement . 

4. The optimal fractional coverage for the wire- wrapping of 
MHT Korodense tubes is -=5% (P w /D w «=20) when used in a 
single tube arrangement. 

5. The thermal conductivity of the tube material effects 
the enhancement characteristics of wire-wrapped tubes. 

6. Some relationship appears to exist between fractional 
coverage and the inside heat transfer coefficient. 

7. The modified Wilson plot data reduction technique 
compensates for vapor velocity effects when using the 
Nusselt theory in determining the outside heat transfer 
coefficient . 

B. RECOMMENDATIONS 

1. Manufacture and test additional wire-wrapped LPD 
Korodense tubes with smaller wire diameters ranging from 
0.1 mm to 0.4 mm to validate the trends found for the 
MHT Korodense tube. 

2. Manufacture and test additional constant pitch wire- 
wrapped copper smooth tubes with wire diameters greater 
than 1.6 mm and less than 0.5 mm and determine the 
optimal wire diameter for specific wire pitches. 

3. Test a set of wire-wrapped corrugated copper or 
copper/nickel tubes to determine the possible effect of 
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tube material on condensation, when compared to the 
titanium corrugated tubes. 

4. Develop a numerical model to approximate the outside 
heat transfer characteristics for corrugated and wire- 
wrapped corrugated tubes. 

5. Fabricate a set of corrugated tubes with outside 
electrical heater elements and conduct tests to validate 
the accuracy of the inside heat transfer correlations 
available for processing data on corrugated tubes with 
and without a HEATEX insert. 

6. Input the Wolverine inside heat transfer correlation to 
the existing program. Reprocess selected corrugation 
tube data and compare results with those obtained by the 
Petukhov- Popov correlation. 

7. Test wire-wrapped corrugated tubes in a tube bundle 
arrangement. Determine if additional overall bundle 
enhancement is possible with inundation effects 
included, and if the optimal fractional coverage shifts 
from that found for single tube tests. 

8. Fabricate and test a set of multistart corrugated tubes 
with and without wire- wrapping to determine if 
additional enhancements can be realized using the 
multistart corrugation geometry. 

9. Conduct additional tests on wire-wrapped smooth and 
corrugated tubes to determine if a relationship exists 
between fractional coverage and the inside heat transfer 
coefficient . 

10. Test the enhancement effects of various wire materials 
when used for wire -wrapping smooth and corrugated tubes. 

11. Investigate the possibility of bonding the wire-wrap to 
the tube surface and test the enhancement effects for 
smooth and corrugated tubes . 
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APPENDIX A. SYSTEM CORRECTIONS 



A. FRICTIONAL TEMPERATURE CORRECTIONS 

Coolant flowing through a tube experiences a temperature 
increase due to frictional heating. The magnitude of the 
temperature rise depends on the internal tube geometry and the 
velocity of the coolant. With the cooling water sump and the 
test section at ambient temperature, measurements were taken 
to determine the temperature increase due to frictional 
effects at flow rates ranging form 10% to 100% on the 
rotameter. Two separate runs were conducted for the MHT 
Korodense tube, one with the HEATEX insert installed and the 
second with no insert. Figure A. 1 show the results for the 
two runs. The data was curve fitted to a third order 
polynomial as shown in Table A. 1. 
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Figure A.l Frictional Temperature Rise for the MHT 
Korodense Titanium Tube with and without a HEATEX Insert 
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Flow Velocity (m/s) 



Table A. 1 FRICTIONAL TEMPERATURE RISE EQUATIONS FOR MHT 
KORODENSE TUBE 



Inside 

Condition 


Polynomial Equation 


HEATEX 

Insert 


T rise = 2 . 564x10 " 4 + 6.263 x10* 4 V w 

+ 2 . 603x10 " 3 ( V ) 2 + 7 . 83 0x10" 

6 (vj 3 


No 

Insert 


T nse = 6.790X10' 5 + 1.632xlO' 3 V w 

+ 8 . 4 09x10 " 4 ( V ) 2 + 1.111x10' 

4 (V w ) 3 


where: T rise = Temperature Rise (K) 

V w = Fluid Velocity (m/s) 
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APPENDIX B. UNCERTAINTY ANALYSIS 



All derived quantities have an inherent uncertainty- 
associated with them when their results are obtained from 
measured values. The uncertainty arises due to inaccuracies 
in the measurement of the physical quantities used in their 
calculations. The difference between the actual value and the 
measured value depends on the accuracy of the measuring 
device, the calibration of the device, and the experience of 
the operator. Eventhough the error associated with any single 
measurement may seem insignificant, it may have a profound 
effect on the accuracy of the derived quantity due to the 
calculations made. 

Kline and McClintock [Ref. 67] developed a method to 
estimate the uncertainty in the derived quantity, R, which is 
a function of several measured variables, x x , x 2 , ..., x^, as 
follows : 



w 




+ 



' dR 
, dX 2 




+ 



+ 





(B.l) 



where : 

W R is the uncertainty of the desired dependent variable 

x x , x 2 , ..., x^ are the measured independent variables 

W 1; W 2 , ..., W n are the uncertainties in the measured 

variables 
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Georgiadis [Ref. 68] presents a complete discussion on the 
uncertainty analysis used for this experiment. Based on 
Georgiadis' developments, Mitrou [Ref. 14] wrote a program for 
calculating the uncertainties in the mathematically derived 
quantities obtained from data collected using the apparatus 
utilized in this research. Modifications were made to the 
program which allow for the data collected on MHT Korodense 
tubes to be processed. Sample uncertainty results for the 
following data collection conditions are presented in this 
appendix. 

1. Vacuum pressure with no insert and high flow rate (3.83 
m/s) . 

2. Vacuum pressure with no insert and low flow rate (1.03 
m/s) . 

3. Vacuum pressure with HEATEX insert and high flow rate 
(3.83 m/s) . 

4. Vacuum pressure with HEATEX insert and low flow rate 
(1.02 m/s) . 

5. Atmospheric pressure with no insert and high flow rate 
(3 . 82 m/s) . 

6. Atmospheric pressure with no insert and low flow rate 
(1.02 m/s) . 

7. Atmospheric pressure with HEATEX insert and high flow 
rate (3.82 m/s) . 

8. Atmospheric pressure with HEATEX insert and low flow rate 
(1.02 m/s) . 
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1. DATA FOR THE UNCERTAINTY ANALYSIS: 



File Name: LMNTNV8 

Pressure Condition: Vacuum 

Vapor Temperature 
Ulater Flow Pate ^ % ) 

Later Velocity 
Meat Flu/. 

Tube-metal thermal conauc. 
Petukov -Popov constant 



48.574 

80.00 

3.83 

1 . 48 1 E+05 
31.0 
3.3857 



Deg C > 



m/ s ) 

U / m 3 ) 

( U/m . K ) 



UNCERTAINTY ANALYSIS : 
VARIABLE 

Mass F Low Rate, Md 
Reynolds Number, Re 
Heat Flu/, q 
Log-Mean- r em Diff, LMTD 
Wail Resist ance , Rw 
Overa II H . T . C . , Uo 
Water-Side H.T.C. , Hi 
Vapor-Side H.T.C. , Ho 



PERCENT UNCERTAINTY 

0.80 
I .09 
1 .87 
1 .63 
4.97 
2.48 
7.47 
6.77 



2. DATA FOR THE UNCERTAINTY ANALYSIS: 

File Name: LMNTNV8 

Pressure Condition: Vacuum 



Vapor Temperature 


= 48.691 


( Deg 


Water F low Rate ( % ) 


= 20.00 




Water Velocity 


= 1 .03 


( m/s ) 


Heat Flux 


= I.149E+0S 


( W/m*2 


Tube-metal thermal conduc . 


= 21.0 


( W/m . 


Petukov- Popov constant 


= 2.3857 





UNCERTAINTY ANALYSIS: 
VARIABLE 

Mass F low Rate , Md 
Reynolds Number, Re 
Heat Flux, q 
Log-Mean-Tem Diff, LMTD 
Wall Res l stance , Rui 
Overall H.T.C. , Uo 
Water-Side H.T.C. , Hi 
Vapor-Side H.T.C., Ho 



PERCENT UNCERTAINTY 

2.99 

3.08 

3.07 

.56 

4.97 

3.13 

7.81 

13.98 
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3. DATA FOR THE UNCERTAINTY ANALYSIS: 



File 


Mane : 


LMNTHU 8 






F^es 


sure Condition: 


1 J a c u u n 






Uapc 


r r enpe r ature 


= 


48.518 


' Deg 


Late 


r FIouj Rate ( X \ 




80.00 




W a t e 


r Ue i oc i t y 


= 


3.83 


( m / 5 


H ea t 


c lu- 




1 .S26E + 05 


U/m 2 


T uOe 


- n e t a 1 trernal 


conduc . = 


21 .0 


' U / m . 



Petuf ov-Popov constant = 3.0112 



UNCERTAINTY ANALYSIS: 
VARIABLE 

Mass R 1 c w Rate, M d 
Reynolds Nunoer , Re 
Heat Flu.-, q 
Log-Mean-Ten 0 i f f , LMTD 
W a 1 1 Resistance, R w 
Overall H.T.C., Uo 
Wa t er-5 i de H.T.C. , Hi 
Uapor-Side H.T.C. , Ho 



PERCENT UNCERTAINTY 

0.80 
t .09 
1 . 83 
1 .58 
4.97 
2.42 
5.92 
6.35 



4 . DATA FOR THE UNCERTAINTY ANALYSIS: 



File Name: 

Pressure Condition: 
U a p o r Temperature 


LMNTHU8 

Uacuun 


48.551 


( Oeg C ) 


Water Flow Rate ( % ) 


= 


20.00 




Water Uelocity 




1 .02 


( m / s ) 


Heat Flu.- 


= 


1 . 202E+05 


( U/m '2 ) 


Tube-metal thermal 


conduc . 


21.0 


( W / m . K ) 


P e t u t ov-Popov const 


art 


= 3.0112 





UNCERTAINTY ANALYSIS: 

VARIABLE PERCENT UNCERTAINTY 



Mass Flow Rate, Md 


3.00 


Reynolds Nunoer, Re 


3.09 


Heat. Flu/. , q 


3.08 


Log-Mean-Ten Diff, LMTD 


.54 


Wall Resistance, Rw 


4.97 


Overa 11 H.T.C., Uo 


3.12 


Water-Side H.T.C., Hi 


6. 36 


Uapor-Side H.T.C., Ho 


11.54 
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5. DAT A FOR THE UNCERTAINTY ANALYSIS: 



Fi le Name: LMNTNA8 

Pressure Condition: Atmospheric 

Vapor Temperature - 

'Jater Flow Pate ( y . ) = 

Water Velocity' 

Heat FIua 

Tube-metal thermal ccrduc . 

P e t u o v - P o p o v constant 



101 kPa > 
99.954 
80 . 00 
o . d ^ 

3. 359E+05 

:i .0 

2.4659 






^ ' s ) 
Urn 2 
1 W / m . r. 



UNCERTAINTY ANALYSIS: 
UARIAELE 

Mass Flow Pate, M d 
Fey nc Ids Number , Re 
Heat F 1 u /* , q 
Log-Mean - r em (Jiff, L M T D 
Wall Resistance, P w 
Overall H . T . C . , Uo 
Water-Side H.T.C. , Hi 
Uapor-Side H.T.C., Ho 



PERCENT UNCERTAINTY 

0.S0 

1.11 

1.10 

.61 

4.9^ 

1 .26 
7.69 
4.97 



6. OATA FOR THE UNCERTAINTY ANALYSIS: 



File Name: LMNTNA3 

Pressure Condition: Atmospheric 

Vapor Temperature 
Water Flow Rate (*£) 

Water Velocity = 

Heat Flux 

Tube-metal thermal conduc . = 

Pe t uk o v- Pop O', constant 



(101 k P a ) 
39.981 
20.00 
1 .02 

3. I I4E + 05 
21 .0 
2.4659 



( Deg C ) 

( m / s ; 

( W/m A 2 ) 

( W / m . K ) 



UNCERTAINTY ANALYSIS: 
VARIABLE 

Mass Flow Rate , Md 
Reynolds Number, Re 
Heat Flux, q 
Log-Mean-Tem Diff, LMTD 
Wall Resistance, Rw 
Over a 1 1 H.T.C., Uo 
Water-Side H.T.C. , Hi 
Vapor-Side H.T.C. , Ho 



PERCENT UNCERTAINTY 

3 . 00 
3.10 

3.04 

.21 

4.97 

3.05 
8.04 

13.03 
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7. DATr C 0F T-E jNCERT^INTV FNAL>SIS: 



File Mane: LMNTHAS 

Fressjre Condition: Atmospheric 

Vapor Temperature = 

Wate r Flow Rate ‘ ’i ) 

Wa t er Ue i oc 1 t y = 

Heat F i u a 

Tuoe-metal f h e r m a 1 conduc. - 
P e t u ► ov -Popov constant 



M 0 1 r Pa 
93.924 
80 . 00 
3.32 

4.083E+05 
21 .0 
3.3112 



Deg C 



' m / S ; 

W/m 2 
( W / m . K 



UNCERTAINTY ANALYSIS: 
VARIABLE 

Mass Flow Rate, Md 
Reynolds Number, Re 
Hea t F I u> , q 
Log-Mean-Tem Diff, LMTD 
Wa 1 1 p es i s t a nee , Rw 
Over a 1 1 H . T . C . , Uo 
Water-Side H.T.C. , Hi 
Vapor-Side H.T.C. , Ho 



PERCENT UNCERTAINTY 

0.80 
1.11 
1 .09 
.59 
4 . 97 
1 .24 
5.68 
4.25 



8. DATA FOR THE UNCERTAINTY ANALYSIS: 



File Name: LMNTHA8 

Pressure Condition: Atmospneric 

Vapor Temperature = 

Water Flow Rate (7.) = 

Wa t er Ve l oc i t y = 

Heat Flu-. = 

Tube-metai thermal conduc. = 

Pet ut ov-Pcpov constant = 



(101 k Pa ) 
99.999 
20.00 
1 .02 

3 . 367E + 05 
21 .0 
3.31 12 



( Deg C ) 

( m / s ) 

( W/m '2 ) 

( W/m.K ) 



UNCERTAINTY ANALYSIS: 
VARIABLE 

Ma S5 Flow Rate, Md 
Reynolds Number, Re 
Hea t Flux, q 
Log-Mean-Tem Diff, LMTD 
Wall Res i 5 tance , Rw 
Over a 11 H.T.C., Uo 
Water-Side H.T.C. , Hi 
Vapor-Side H.T.C., Ho 



PERCENT UNCERTAINTY 

3.00 
3. M 

3.04 
.19 

4.97 

3.05 
G. 14 
9.73 
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APPENDIX C. DRPML PROGRAM LISTING 

The data collection and reduction program, DRPML, which 
was used throughout this study follows: 
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1 000 I 
1 00b 1 
1 007 I 
1903 » 
1003' 
10101 
101 Z I 
10131 
1014 1 
1015' 
1017' 
10181 
1013» 
1020 ' 
1 02 1 ' 
1022 ' 
1023' 
I 024 f 
1025' 
1026 i 
10271 
1023* 
1029' 
1030’ 
1031 i 
1 032 i 

I 8 j j i 

I 834 ' 
1835 
1036 
I 037 
1 033 
1039 



CRPML LONG i 

REVISED FFQM DRPC- : SEP 990 0 ’ c.EEFE / 

r j BE USED WITH NuN- INSTRUMENTED TUBES ONL < 

T~c.ES OnTA IN THE FORMAT QF S'JENSEN/O ’ KEEFE- LONG 
CAM FEPFOCESS ANf MON- INSTRUMENTED DATA 

THIS PROGRAM WAS USED TO COLLECT ALL THE NON- 

INSTFUMENTED DATA TAKEN BY L0N6 ( JAN-JUN 1993) FOR TITANIUM TUBES 

MEANING OF ALL FLAGS IN PROGRAM 

IFT : FLUID TY P E 

ISO: OPTION WITHIN PROGRAM 

IM: INPUT MODE 

IUIL: VALUE OF Ci USED 

IFG: FINNED OP SMOOTH 

INN: INSERT T /PE 

I U T : LOOP NO. WITHIN PROGRAM 

I UTH : ALTERNATIVE CONDENSER TUEES 

iric : TUBE MATERIAL 

I TDS : TUBE DIAMETER 

IFC : FRESSURE CONDITION 

INF: DIMENSIONLESS FILE REQUIRED 

IFF: FLOT FILE REQUIRED 

I OV : OUTPUT REQUIRED 

I H I : INSIDE HTC CORRELATION 

IOC: OUTSIOE HTC THEORY /CORRELATION 

COM /Cc/ C ( 7 ) 

COM /CcE5/ T5S( 5 ) 

COM /CcSS / TS6 ( 5 ) 

COM / Ccb7 / T57 ( 5 ) 

COM ■ CcS3/ T58<5! 



1040 COM 'Fid/ Ift.Istu 
104 1 DIM Enf ( C0 ) ,Tw( G ) 

1042 COM -Pr/ Qpa ( 42 ) , T f m( 42 ) , T f nr , I pc , Qpr 

1043 .COM / W 1 1 / Nrun , 1 1 m , Iwt h , Inc , I fe , I j ob , I wd , I f g , Ipco , I f t o , I w i l , I h l , I oc , I nan ,K 



cu ,Re.- p ,Rm , A.-. 



1 044 
1045 
1047 
1 048 
1 049 
1 050 

1051 

1052 



COM /Geon/ 01 ,D2 .Dx ,Do ,L , L 1 ,L2 

COM /Fnc/ I s t uo , I nn , I w t ho , Vw , Vus , I nano 

DATA 0.10086091 ,25727.9436 9,-767345.8295,78025595.81 

DATA -9247486 539 ,6. 97683E+1 1 ,-2 . 68 1 92E+ 1 3 ,3 . 94078E+ 1 4 

READ C< * ) 

DATA 273. I5.2.5943E-2 , -7 . 267 1 E-7 ,3 . 294 1 E- 1 1 ,-9 . 77 1 9E- 1 6 ,9 . 7 1 2 1 E-20 
READ T55( * ) 

DATA 273. 15 ,2 .5978E-2 ,-S. 9853E-7 ,-3. 1242E-1 I ,1 .3275E-1 4 ,-1 .0188E-1 8 
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1053 

1054 

1055 

1056 

1057 

1 058 

1055 

1050 

1 06 I 

1062 

1063 

1 064 

1065 

1066 

1 067 

1069 

1090 

1 106 

1 1 08 

Mil 

1112 

1115 

1116 

1117 

1118 

1119 

1 1 20 

1 123 

1 1 24 

1 126 

1 129 

1 132 

1 133 

1 1 35 

1 136 

1 1 38 

1141 

1 1 44 

1 1 45 

1 1 46 

1 1 48 

1 150 

115 1 

1 153 

1 156 

1171 



READ T56 ( * ) 

DATA 273.15 .2.5923E-2 ,-7.3933E-7 ,2.8625E-1 1 
READ T57( ♦ ) 

DATA 273. 15 ,2. 5 93 IE-2 ,-7.5232E-7 .4.0657E-1 1 
P E AD T 5 8 ( * ) 



,1 .9717E-15.-2.2486E-19 



,-1 .2791E-15 ,6 . 4402E-20 



Dr=. 015375 1 Outside diameter of the outlet end 

Ds5p=.1524 ' Inside diameter of stainless steel test section 



Aa=PI*Dssp 2 / 4 



A 1 p 2 = 0 . 

L= . 13335 
L1 = . 060325 



Condensing length 
Inlet end "fin length 



L2 = . 034925 1 Outlet end "fin length" 

PRINTER IS 1 

BEEP 

PRINT USING " 4'X , ' "Select option: 

PRINT USING "6X."" 0 Take data or re-process 

PRINT USING " 6 X , " " 2 WILSON Analysys 

INPUT Iso 

I SG= I 501- 1 

IF Iso 1 THEN 3094 



previous 



da t a 



EEEP 

INPUT "SELECT FLUID C 0=WATEP , 1=R-113, 2 = EG) M ,Ift 
I f t o= I ft 
EEEP 
I j ob = 0 

INPUT "ENTER INPUT NODE ( 0 = 3054A , 1 =F I LE ) " , I m 

I m= i m+ 1 

BEEP 

IF I m= 1 THEN 

INPUT "ENTER MONTH, DATE AND TINE t MM : DO : HH : MM : S3 ) M .Date* 
OUTPUT 709; "TO" ;Oate$ 

OUTPUT 709; "TO" 

ENTER 709;Date$ 

END IF 

IF I j od = 1 THEN 
BEEP 

INPUT "SKIP PAGE AND HIT ENTER “ ,0k 
END IF 

PRINTER IS 701 
IF I rn= 1 THEN 

ENTER ^09;Oate$ 

PRINT " Month, date and time :";Date$ 

END IF 
PRINT 

PRINT USING "10X I ""NOTE: Program name : ORPML 

IF I job = l THEN 1 189 
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1 74 


BEEP 






1 96 


INPUT SELECT < C l : 0 = F I NO , 1 “STORED Ci/',Iwl 




1 39 


IF ]>=1 THEN 






i b 3 


BEEP 






i 35 


INPUT " G I 


VE A NAME FOR THE RAW DATA FILE" ,0_f 1 1 e$ 


1 98 


PRINT USING '16* .""File name : " " , 1 J A “ ; D_f l 


leS 


301 


CREATE BDAT D_f l Ie$ ,30 




204 


ASS I SN @F 


lie TO G_ f 1 1 e $ 




207 


BEEP 






Z 1 0 


INPUT “ENTER GEOMETRY CODE ( 1 =FINNED ,0 = PLAIN )“ 


, I fg 


21 1 


I nn = 0 






Z 1 z 


PRINTER I 


S 1 




Z 1 6 


BEEP 






Z 1 7 


PRINT " 


ENTER INSERT TYPE: ' 




218 


PRINT " 


0 = NONE < OEFnULT ) ” 




Z 1 9 


PRINT " 


1 = T i J I S T E D TAPE" 




ZZ0 


PRINT “ 


2 = U I RE URPiP" 




ZZ 1 


PRINT " 


Z=HEPTEX 11 




ZZ2 


INPUT Inn 






2Z5 


OUTPUT OF 


ile;Ifg,Inn 




n ^ n 


I w t =0 1 


FOR UN I NSTRUNENTEO TUBE 




ZZ8 


Fh=0 






23 1 


Fp = 0 






Z34 


F ui- 0 






Z35 


I 5 t u = 0 






236 


I 5tUO= I S tu 




238 


IF I fg=0 


THEN 1241 




239 


INPUT "FIN PITCH, HEIGHT AND WIDTH IN MM, Fp ,Fh ,Fw 


24 1 


OUTPUT ®F 


ile;Iwt,Fp f Fui,Fh 




242 


ELSE 






249 


BEEP 






250 


PRINTER IS 1 




252 


PRINT " 


STUDENT’S DATA TO BE REPROCESSED:" 




LD 
f J 


PRINT ' 


0 = SWENSEN/O ’KEEFE/LONG ( DEFAULT )" 




254 


PRINT " 


1 =UAN PETTEN/'MITROU/COUMES/GUTTENOORF 


" 


255 


INPUT 1st 


u 




25S 


I 5 t UO= I 5 t U 




253 


BEEP 






259 


PRINT 






260 


IF I stu= 1 


THEN 




26 1 


PRINT 


" STUDENT NAME : " 




262 


PRINT 


“ 0=UAN PETTEN" 




263 


PRINT 


” 1 =M I TROU " 




264 


PRINT 


" 2=COUMES " 




265 


PRINT 


” 3=GUTTEND0RF " 




266 


ELSE 
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1 2B7 PRINT •' 4 = StJENSEN" 

1258 PRINT " 5=0 ’ KEEFE " 

1269 PRINT •' 6 =LGNG" 

1271 END IF 

1272 INPUT I nan 

1273 I nano = I nan 

1275 BEEP 

1275 INPUT " 6 1 UE THE NAME OF THE EXISTING DATA F IlE" ,D_f i lei 

1277 PRINTER IS 701 



1278 
EN 


IF 


I nan=0 


THEN 


FRINT 


USING 


“ 1 6X . 


/' "Data 


t ak en 


by 


: VAN PETT 


1 279 


IF 


I nam = 1 


THEN 


PRINT 


USING 


” I E X , 


"Data 


taken 


by 


: MITROU"" 


1 280 


IF 


I nam = 2 


THEN 


PRINT 


USING 


M 1 6 X 3 


.""Data 


taken 


by 


: COUMES"" 


1281 
RF 


IF 


I nam-3 


THEN 


PRINT 


USING 


” 1 GX , 


.""Data 


taken 


by 


: GUT7END0 


1282 


IF 


I nam = 4 


THEN 


F INT 


USING 


1 6 X , 


.""Data 


t ak en 


by 


: SWENSEN" 


1283 


IF 


I nam=5 


THEN 


PRINT 


USING 


” 1 6 X a 


, " "Data 


taken 


by 


: O’KEEFE" 


1284 


IF 


I nam = 6 


THEN 


PRINT 


USING 


M 1 6 X , 


.""Data 


taken 


by 


: LONG 



1285 PRINT USING "!6X,""This analysis done on file : , I 0A " ; D_f l 1 e$ 

1 2S6 PRINTER IS 1 

1287 BEEP 

1288 INPUT "ENTER NUMBER OF DATA SETS STORED" ,Nrun 

1289 ASSIGN QFile TO D_file$ 

1290 ENTER @F i le ; I f g , Inn 

1291 IF I s t u = 0 THEN 

1292 ENTER @F i le ; Iwt ,Fp ,Fw ,Fh 

1293 ELSE 

1294 IF I f g=0 THEN ENTER SFileslut 

1295 IF I f g= 1 THEN ENTER @Fi 1 e ; Fp ,Fw ,Fh 

1 296 END IF 

1297 END IF 

1298 IF I j ob= 1 THEN 1537 

1299 IF Ift .0 THEN 1345 

1300 BEEP 

1301 PRINTER IS 1 



1302 


PRINT 


USING 


" 4X , 


/'"Select tube t ype : 






1303 


PRINT 


USING 


"GX. 


/' "0 Copper Tube' 






1 30G 


PRINT 


USING 


"GX 


/‘"I Wolverine Korodense LPD 


T itanium 


Tube 


1 320 


PRINT 


USING 


"GX . 


/' " 2 Smooth Titanium Tube""" 






1 321 


PRINT 


USING 


" 6X 


, " " 3 Wolverine Korodense MHT 


Titanium 


Tube " " " 


1323 


INPUT 


I wt h 










1324 


I w t ho : 


= I ut h 
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1526 

1323 

1329 

1 330 

I 33! 

1332 

1 333 

1 334 

1 335 

1 336 

1337 

1 336 

1 u 0 3 

1 340 

! 34 1 

1 342 

1 343 

1345 

1343 

1349 

1351 

1352 

1357 

1360 

1 363 

1366 

1367 

1369 

1 372 

1375 

1373 

1383 

1384 

1 399 

1402 

1 405 

1457 

1458 

1 459 

1460 

1461 

1462 

1463 

1 464 

1 466 



BEEP 










PRINT 


USING 


■U 




elect tube Enhancement used:" 


PRINT 


USING 


' b ' , 


0 


SMOOTH TUBE 


rS I N T 


USING 


' 6 .< , 


' ! 


PINNED TUBE 


PRINT 


USING 


"6a , 


• 2 


WIPE-WFAPPED SMOOTH TUBE 


FPINT 


USING 


' 6 < , 


« ~7 


LPD h OPOCENSE TUBE ' ' " 


PRINT 


USING 


" 6 X t 


"4 


WIFE-WRAPPED LPD KORODENSE TUBE"" 


PR INT 


USING 


"6X , 


" 5 


MHT KOPODENSE TUBE 


PRINT 


USING 


" 6X 


1 "6 


WIRE-WRAPPED MHT KORODENSE TUBE"" 


INPUT 


I typ 








PRINTER 15 ' 


701 






BEEP 










PRINTER IS 


1 






PRINT 


USING 


■’ 4X , ' 


'"Select Material Code '• " " " 


PRINT 


USING 


" 6 X 


• - a 


Copper 1 Stainless steel" 11 ’ 


PRINT 


USING 


"6a,' 




Aluminum 3 90:10 Cu-Ni" " " 


PRINT 


USING 


“ 6X , 


‘ " 4 


Titanium " “ " 


INPUT 


Inc 








PRINTER IS 


1 






BEEP 










1 1 ds= : 


1 








IF Iut h=0 THEN 







PRINT 


USING 


" 4X , 


.""SELECT TUBE DIA TYPE: 


PRINT 


USING 


" 6X , 


, 0 


SMALL 


PRINT 


USING 


" 6X , 


, “ 1 


MEDIUM (DEFAULT ) 


PRINT 


USING 


" 6X , 


.. .. ^ 


LARGE 


INPUT 


I ids 









END IF 

PRINTER IS 701 
IF I w t h = 0 THEN 

D i = . 0 1 27 i ID OF MEDIUM 8, LARGE TUBES 

Do=. 01905 * OD OF MEDIUM TUBE 

END IF 

IF I uj t h = 1 THEN 
Di=.01 347 
Do= .01 585 
END IF 

IF Iut h = 2 THEN 
Do= . 0 1 595 
Di= . 0 1 386 
END IF 

IF Iut h=3 THEN 
Do=.01 607 
Di=. 01353 
END IF 
D1 =.01 905 
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I 467 

I 468 

I 469 

1 471 

I 474 

t 477 

I 478 

1 479 

1 484 

I 485 

1 487 

1 489 

I 492 

1 495 

1 497 

I 498 

1501 

1504 

1507 

' 508 

1510 

1537 

1538 

1 543 

I 546 

1549 

1552 

I 556 

1557 

1559 

1559 

1560 

1561 

1562 

1563 

1564 

1565 

1567 

1568 

1569 

1570 

1 571 

1572 

1 573 

1574 

I 575 









D2=. 01537 
IF ltds=0 THEN 
Do= .0127 
Di = . 003525 
END IF 

IF 1 1 ds = 2 THEN Do=.025 

IF I w t h= 1 OR I w t h=3 THEN D 1 ■= . 01535 

IF Iwth=1 OR Iuth=3 THEN D2=. 01535 

IF Iwth=2 THEN Di=. 01537 

IF Iwt h=2 THEN 02= . 01587 

IF I me =0 THEN Kcu=335 

IF I m c = 1 THEN K c u = 1 6 

IF I mc = 2 THEN Kcu=l67 

IF Imc = 3 THEN Kcu=45 

IF I mc = 4 THEN Kcu = 2 1 .0 

Pm = Do*L0G(Do/DiJ/(2*Kcu W a 1 1 resistance based on outside area 

BEEP 

INPUT "ENTER PRESSURE CONDITION 0 = U , 1 =A ) " f I pc 
I p c o = Ipc 
Inf =0 



BEEP 



I f e = 1 



PRINTER 13 701 

PRINT USING " 1 5 X , " " Th i s analysis includes end-fin effect 
PRINT USING " 1 6X , " " Therma 1 conductivity = " M ,3D.D," 
PRINT USING " 1 6 X , " " I ns l de diameter, Di = n ",DD.DD, 
PRINT USING " 1 BX , " " Ou t s i de diameter, Do = " " ,DD.DD, 
I h i =0 



( W/m .K ) 
( mm ) " " 
( mm ) " " 



PRINTER IS 
PRINT •’ 
PRINT " 
PRINT " 
PRINT " 
INPUT Ihi 



SELECT INSIDE CORRELATION:" 
0 = 5 IEOER-TATE ( DEFAULT ) 
1 =5LE I CHER -ROUSE " 

2 =PETUK HOU-P'OPOU " 



IF Ih l =0 THEN 
BEEP 

INPUT " SELECT REYNOLDS EXPONENT 
END IF 



,Rexp 



I oc=0 
EEEP 
PRINT 
PRINT " 
PRINT " 
PRINT " 
INPUT Ioc 
EEEP 



SELECT OUTSIDE THEORY/CORRELATION FOR WILSON ANALYSIS: 
0=NUSSELT THEORY (DEFAULT)" 

1 =FU J 1 1 (1979) CORRELATION" 












" “ ; Kcu 
; D l ♦ 1 000 
; Do* 1 000 
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1576 It ii=1 

1577 PRINT 

1575 PRINT SELECT COOLANT TEMPERATURE RISE MEASUREMENT: 

1579 IF I 5 t u ”0 THEN PRINT " 0=SIN6LE TEFLON T/C" 

1550 PRINT " 1 =OUAFTZ THE c MOMETE° ' OEFAULT i " 

1551 D P I NT 2= I 0-JUNCT ION THERMOPILE" 

1 58Z INPUT I tm 

1593 PRINTER IS 701 



1 594 


IF 1 1 n=0 


THEN PRINT 


USING " 


1 6X 


,""This analysis uses the SINGLE "EFLGN T C 


readings ' " " 










1555 


IF 1 1 m = 1 


THEN PRINT 


USING " 


I6X 


/'"This analysis uses the QUARTZ THERMOMETER 


readings" " " 










1556 


IF 1 1 n = 2 


THEN PRINT 


USING " 


16X 


/‘"This analysis uses the 10-JUNCTION THERMO 


PILE 


reaGings " 


" " 








1537 


I 1 c - 1 1 


FOP MODIFIED WILSON 






1 553 


IF I h i =0 


THEN C l = . 0 








1555 


IF I h i = 1 


THEN C l = 1 . 








1596 


IF I n i =2 


THEN C ; = 1 . 








1559 


IF I w i 1 = 


1 THEN 








1 601 


BEEP 










1 602 


INPUT 


"ENTER Ci 


IF DIFFERENT FROM STORED UALUE",Ci 


1603 


END IF 










1 604 


PRINTER 


IS 701 








1 605 


IF I hi =0 


THEN PRINT 


USING " 


1 6X 


/'"Modified Sieder-Tate coefficient = "",Z 


.40" 


; C i 










1606 


IF I h i =0 


THEN PRINT 


USING “ 


1 6X 


/'"Chosen Reynolds No. exponent = ,G 


.DO" 


; R e a p 










1607 


IF I h l = 1 


THEN PRINT 


USING " 


1 6X 


/‘"Modified S 1 e i c her-Rous e coefficient 


" .z.< 


40" ; C l 










1609 


IF I h i-2 


THEN PRINT 


USING " 


1 6X 


/‘"Modified Pe t uk hov-Popov coefficient 


" . Z . ■ 


4D " ; C l 










1609 


IF I nn=0 


THEN FRINT 


USING " 


1 6X 


/'"Using no insert inside tube 


1610 


IF I nn=2 


THEN FRINT 


USING " 


1 6X 


/‘"Using wire wrap insert inside tube 


1611 


IF Inn=3 


THEN FRINT 


USING “ 


1 6X 


/‘"Using HEATEX insert inside tube 


1612 


IF I s t u = 1 


0 THEN 








1613 


IF I nn= 1 THEN PRINT USING " 


16X /'"Using twisted tape insert inside tube" 


1614 


ELSE 










1615 


IF I nn= 1 THEN PRINT USING " 


16X /'"Using wire wrap insert inside tube 


1616 


END IF 










1619 


IF I ic=0 


AND I f e= 1 


THEN Ac= 


26. 


4 


1630 


IF I ic- 1 


THEN Ac=0 . 








1621 


BEEP 










1622 


IF I j cb = 


1 THEN 1643 








1623 


PRINTER 


IS 1 









1624 1 INPUT "NAME FOR A TEMPORARY PLOT FILE (TO BE PURGED )" ,P_f 1 1 e$ 
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1625 P_f l leS="DUMMY" 

1626 SEEP 

1634 CREATE BDA T R_ f ils'£,10 

1644 DESIGN ?Fiiep TO P_ f ile$ 

1 6 4 3 it~ Ijod = ' then 

1651 Io - = 1 

1654 GOTO 1667 

1657 END I c 

1660 BEEP 



1 661 


INPUT "SELECT 


OUTPUT 


1 0=5HORT , 1 


= LONG > 1 


, I 0 V 


1666 


I ov 


= I 0 '/ + 1 














1667 


PRINTER I 


S 701 












1 672 


IF 


I t yP=0 


THEN 


P*INT 


USING 


” 1 6 X , 


"Tube 


Enhancement 


1 673 


IF 


I t / p = 1 


THEN 


PRINT 


USING 


" 1 6X , 


, " "Tube 


Enhancement 


1 674 


IF 


1 1 y p = 2 


THEN 


PRINT 


USING 


" 1 6X , 


, " " T ube 


Enhancemen t 


TUBE" 


















1675 


IF 


i typ = 3 


THEN 


PRINT 


USING 


" 16X ; 


, " "Tube 


Enhancement 


1675 


IF 


Ityp-4 


THEN 


PRINT 


USING 


" 1 6X , 


, " "Tube 


Enhancement 


RODENSE TUBE 














1675 


IF 


1 1 yp = 5 


T HEN 


PRINT 


USING 


" 16X , 


, " " Tube 


Enhancement 


1660 


IF 


I typ=6 


THEN 


PRINT 


USING 


" 16X , 


"Tube 


Enhancement 



PODENSE TUBE 
1632 BEEP 



1 633 


IF 


I m c = 0 


THEN 


P° I NT 


USING 


" 1 6X 


, " "Tube 


material 


1684 


IF 


1 me = 1 


THEN 


PRINT 


USING 


" 1 6X ( 


, " " Tube 


ma t er l a 1 


1665 


IF 


I mc=2 


THEN 


PRINT 


USING 


- 1 6 X 


, " "Tube 


material 


1686 


I F 


I me = 3 


THEN 


PRINT 


USING 


"I6X, 


, " "Tube 


material 


1567 


IF 


I me - 4 


THEN 


PRINT 


USING 


" 16X , 


, " " Tube 


ma t er l a 1 


1688 


IF 


I pc =0 


THEN 


PRINT 


USING 


" 16X , 


, " "Pressure condition 


1689 


IF 


Ipc=1 


THEN 


PRINT 


USING 


" 16X , 


.""Pressure condition 



1S90 1 FPINT USING "t6X,""Fin pitch, width, and height inn) 
2.0D" ; Fp , Fw ,Fh 

1691 IF ( Iw i 1 =0 OR Iw 1 1 =2 ) AND Im=2 THEN 

1692 I j ob= I 

I 693 Iwd= I 

1696 CALL W i 1 son( C i ) 

1699 END IF 
1702 J=0 

1712 IF Iov= 1 THEN 
1722 PRINT 



1723 
1 724 


IF I hi= 1 
PRINT 


THEN 

USING 


" 10X / 


‘"Data Uuj 


Uo 


Ho 


Ts 

1725 


Re,<p " " " 
PRINT 


USING 


" 1 0X , 


’ " 4 ( m/ 5 ) 


( W/m 2-K ) 


( U/m 2 



SMOOTH TUBE 

FINNED TUBE 

HIRE-WRAPPED SMOOTH 

: LPD KQPODENSE TijEE" 

: WIRE-WRAPPED LPO KO 

: MHT KORODENSE T| JBE “ 

: WIRE-WRAPPED MHT KO 



COPPER 

STAINLESS-STEEL 

ALUMINUM 

90/10 CU/NI 

TITANIUM 

VACUUM 

ATMOSPHERIC 

“ •' ,DD.DD ,2X ,Z .DD ,2X , 



Op Tc f 

K) ( W/ci A 2 ) <C) 
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I 7 2 5 £ l 5 £ 



”13 


FFINT j-ING 


' o* , 


1 " Da t a 


'vvU 


Uo 


Ho 


Op 


723 


i 5 

C FI NT -SING 


' 1 ax , 


' " & 


"1/ 5 


W / pi 2 - t* 


U) / n 2 ~ t* 


U ' m 



■ C ■ 0 ) 

1730 END IF 
1740 END IF 
1747 Z.*=0 

1750 Z 2 = 0 
1753 Z*.y =0 
1 756 Zy=0 
1759 S,.=0 

1763 S y = 0 
1765 5*5=0 

1768 5-* y = 0 

I 77 1 Go_on= I 
1774 Repeat: 1 
1777 J = J + } 

1780 IF I «i=1 THEN 
1 7S3 BEEP 

1786 INPUT "LIKE TO CHECK NG CONCENTRATION ( 1 =Y ,0=N )7 “ ,Ng 

1739 IF J=1 THEN 

1792 OUTPUT 709 ; " AR AF40 AL4I VR5 " 

1 795 OUTPUT 709 ; "AS SA" 

1798 END IF 

1801 BEEP 

1804 INPUT 'ENTER FLOWMETER READING" ,Fm 

1807 OUTPUT 709; "AR AF50 AL62 VR5 " 

1510 OUTPUT 709; "AS SA" 

1813 ENTER 709;Etp 

1816 OUTPUT 709; "AS SA" 

1319 BEEP 

1.522 INPUT "CONNECT UOLTAGE LINE" ,0k. 

1825 ENTER 739;6vol 

1 828 BEEP 

1531 INPUT "DISCONNECT UOLTAGE LINE" ,0k 

1 334 IF B vo 1 1 . 1 THEN 

1 337 BEEP 

1 640 BEEP 

1843 INPUT "INVALID UOLTAGE, TRY AGAIN i " ,0k 

1346 GOTO 1819 

1849 END IF 

1 858 OUTPUT 709; "AS SA" 

1861 ENTER 709;Bamp 

1862 E t p = Et p * I . E + 6 
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1S63 OUTPUT 709; AF40 AL47 UR5 " 

1374 rjn = 7 

1575 FOP '=0 TO Tin 

1573 OUTPUT T09;"*iS SA " 

1385 5e=0 

I S83 FOR H = I TO 10 

1691 ENTER 709;E 

I 394 Se=5e+E 

IBS' 7 NEXT K 

1900 Enf( I > = A6S(Se/'10 ) 

1916 Enf ( I )=Er>f ( I ) » 1 . E + 5 

1913 NEXT I 

192! OUTPUT 709; "A3 SA" 

1924 OUTPUT 713; " T 1 R2E " 

1927 WAIT 2 

1930 ENTER 7 1 3 ; T 1 1 

1 933 OUTPUT -713; " T2P2E " 

1936 WAIT 2 

1 939 ENTER 7 1 3 ; T2 

1942 OUTPUT 713; "TIR2E" 

1945 WAIT 2 

I 948 ENTER 7 1 3 ; T 1 2 

1951 T 1 =< T I I +T 1 2 ) ♦ . 5 

I 954 OUTPUT 713; " T3RZE " 

I960 BEEP 

1970 INPUT "ENTER PRESSURE 6A6E READING (Pga)",Pga 

1971 Pvap 1 =Pga*6694 . 7 i PSI TO Pa 

1972 OUTPUT 709; "AR AF64 AL64 UR5 " 

1973 OUTPUT 709; "AS SA" i PRESSURE TRANSDUCER 

1974 S s = 0 

1975 FOR K = I TO 20 

1976 ENTER 709;Etran 

1 877 5s = Ss *E t ran 

1 878 NEXT K 

1 879 Ptran=ABS( Ss/20 ) 

1900 BEEF 

1991 i PRESSURE IN Pa FROM TRANSDUCER 

1 987 Pvap2= ( -2 . 93604*Ptran+1 4.7827 ) ♦6894. 7 

1 98E ELSE 

1 98G IF lstu=0 THEN 

1 989 ENTER ®Fi le;Bvoi ,Eamp ,Etp ,Fm ,T1 ,T2 ,Pvap1 , Pvap 2 ,Enf ( * ) 

1990 ELSE 

1992 ENTER @File;Bvoi l Bamp l Utran > Etp I Enf(0),Ernf(1 ),Emf(2) l Enf(3),Enf(4),F 

m , T 1 , T2 ,Phg ,Pwa t er 
1994 END IF 

1996 IF J=1* OR J = 20 OR J=Nrun THEN 
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1997 Mg- 1 

1993 ELSE 

1993 Ng =0 

2000 END IF 

2002 END IF 

2003 IF lstu=0 THEN 

2008 T s t earn 1 =FNT , 5 v57 \ En f > 0 M 

2009 T s t earn 1 =Ts t earn 1 -273 . I 5 

2010 Tsteam2~FNT V 'sv56(EmfM ; ) 

2011 Ts t eam2=Ts t eam2-2*73 . 1 5 

2012 Ts t eam = T 5 1 earn I 

2015 Troom=FNTv5v58(Emf(2M 

2023 Trocm= Trocm- 273 . 1 5 

2038 Tcon=FNTv5/53(Enf'7') 

2033 Tcon=Tcon-2' : 3.15 

2042 ELSE 

204 3 F 5 t e a m = F N T v 5 v ( Em f ' 0 > ) 

2044 Troorn=FNT^svEmf'. 3)) 

204 5 Tcon=FNTvs v ( Ertf < 4 ) ; 

2046 END IF 

2048 Fsat=FNPvst(T steam) 

2050 Rohg= 1 352 3-1 2 2 * ( Tr oom-2 6 .85 )/50 
2053 Rowa t er“FNRhow ( Troon ) 

2063 IF lstu=0 THEN 

2081 Ptest 1 = P v a p 1 

2082 Ptest2=Pvap2 

2083 ELSE 

2084 Ptest2 = ‘. Phg*Rohg-Puater*Ro water '*9.81 / 1 000 

2085 END IF 

2087 PI s=Psa t * 1 .E-3 

2088 PI p = P t es 1 2 * 1 .E-3 

2089 PI* t-Pks 

2092 Tsat = FNTvsp< Psat ) 

2098 Os t = FN0v s 1 1 T s t earn ) 

2104 Ppng=( Ptest 2-Psat ) / P t e 5 t 2 

2121 Pps t= 1 -Fpng 

2122 Nwv= 18.016 

2123 IF I f t = 1 THEN Mwv=137 i TO BE CORRECTED 

2124 IF I f t =2 THEN Mwv=62 

2125 Of ng*< Pt estZ-Psat ) / P t e s 1 2 

2126 M f ng= 1 / ( 1 + ( 1/Ufng-i >*Mwv/29.97> 

2127 M f ng = P1 f ng * 1 00 

2128 Of ng=0f ng* 1 00 
2131 BEEP 

2134 IF Iov=2 THEN 
2137 PRINT 
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2 1 3S * RECORD TIME OF TAKING DATA 

2139 IF In=i THEN 

2140 OUTPUT *709 ; " TD" 

2U1 ENTER 709; T o I d $ 

2142 END IF 

2144 PRINT USING " 1 0X , " " Da t a set numcer 

2145 OUTPUT 709 ; " HR AF40 Al40 UR5 ' 



,DD ,4 < ,1 4 A" ; J , ToldS 



2146 


OUTPUT 709; AS SA" 




2149 


END IF 




2152 


IF Iov=2 AND Ng= 1 THEN 




2 155 


PRINT USING 1 I0X " 


P sa t P t ran 


2158 


PRINT USING “10X,““ (kPa) < KPa ) 


2161 


PRINT USING " 11 X , ]( 3D .DD 


,2/ » ,1 ( 3D . DD , 


PI'P , 


T s t earn , T sa t , M f nq 




2164 


PRINT 




2167 


END IF 




2 1 7 0 


IF M p ng .5 THEN 




2173 


BEEP 




2176 


IF In* 1 AND N g = 1 THEN 




2179 


BEEP 




2182 


PRINT 




2 1 85 


PRINT USING M 10X 


Energire the 


2183 


EEEP 




2 191 


INPUT "OK TO ACCEPT 


THIS RUN < 1 = 


2194 


IF OK =0 THEN 




2197 


BEEP 




2200 


DISP ‘NOTE: THIS 


DATA SET U I L 


2203 


WAIT 5 




2206 


GOTO 1790 




2209 


END IF 




2212 


END IF 




2215 


END IF 




2216 


IF Im= 1 THEN 




222 1 


IF Fm ,10 OR Fn 100 THEN 


2224 


I f m = 0 




n -i -t y 


BEEF 




2230 


INPUT ' INCORRECT FM 


( 1 ^ACCEPT ,0=i 


2233 


IF I f n = 0 THEN 1804 




2236 


END IF 




2239 


END IF 




2242 


i ANALYSIS BEGINS 




2243 


IF lstu=0 THEN 




2252 


Til=FNTvsv58(Emf<3) ) 




2262 


T i2=FNTvsv55( Enf( 5 ) ) 




2272 


Tol=FNTvsv58(Emf( 4 ) ) 




2282 


T o2=FNT vsv5S( Enf < 6 ) ) 





Tmea s 

(C ) 



Ts a t 
(C ) 



NG 

Mo 1 al 



X 

2* >" ;P^ s 



acuum system 
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:3o: 
23 1 : 



2342 



l 1^=1 

To 1 =To1 *273 
T c2=To2-2^3 



'de i 



□ e i 



15 

• 5 
1 c 



= T o 1 -T: 1 
>Tz2-Ti2 



2352 Tde i 3 = T2-T 1 

2353 Etpl=£mf(3)^Etp/20. 

2354 D t de = 2 .5931 E-2 - l . 504S4E-6 *Et p 1 + I . 2 1 70 1 E- 1 0*Etp T'2-5 . i 1 64E- 1 5 



201 E - 


1 3 * E i p r 4 


2355 


Tns=Dtde*Etp/10. 


no 
LD 
m 
t j 


T o3 = T i 1 +Tr i 5 


2359 


IF 107=2 THEN 


2351 


PFINF USING "1 


PILE 





2362 


F D I M T USING ’* I 


2 364 


PRINT USING '2 


2365 


END IF 


2367 


Er 1 =rBS( T i 1 — T 1 ) 


2358 


Er3= A5S ( T i 2 -T 1 > 


2369 


PRINTER IS 1 


2370 


BEEP 



TIM1 TOUT 1 TIN3 T0UT3 DELTI 

■ TEFLON ) t QUARTZ ) 

3D.DD ,2X >' ; Ti 1 ,To1 ,T1 ,T2 ,Tdel 1 , Tde 1 3 ,Tn 



2375 Er 2 = nBS ( v T2 - T 1 ) -( Tr i s ) ) / ( T2-T 1 ) 

2377 IF Er2 .05 AND I 1 THEN 

2378 EEEP 

2379 PRINT "QCT AND T-PILE DIFFER BY MORE THAN 5% 

2380 Qk2= I 

2381 IF 0 k 2 = 0 AND Er2, .05 AND I n= 1 THEN 1730 

2332 END IF 

2383 PRINTER IS 701 

2384 ELSE 



2385 T i =FNT /sv< Em f ( 2 ) ) 

2336 6raa=FNGradMTl+T2 , *.5) 



2387 To = Ti+ABS(Etp)/( 1 0*Grad )* 1 . E+6 

2338 T i 1 = T i 



2339 To3=To 

2391 END IF 

2332 IF lstu=0 AND Itm-0 THEN 

2393 T 1 i = T l 1 

2394 T2o=To 1 



2395 END IF 

2396 IF I t n- 1 THEN 

2397 T 1 i =T 1 

2398 T2q=T2 

2399 END IF 

2400 IF 1 1 m = 2 THEN 



Etp 1 



DELT 
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240 1 

2402 

2403 

2404 

2405 
Z 40b 
2407 
Z408 
2 4 1 0 
241 1 
Z 4 1 Z 
2413 
241 4 
2415 

2417 

2418 
241 g 
2420 

2444 

2445 
2445 

2447 

2448 

2449 

2450 

2451 

2452 

2453 

2454 

2455 
2455 
245 1 
2462 

2465 

2466 

2467 
2458 

2470 

2471 

2472 

2473 

2474 

2475 

2476 

2477 
2476 



T 1 l =T i 1 
T2o = To3 
END IF 

T a v g = ( T 1 i - T2 o / * . 5 
I f t =0 

Cpw = FNCpuM Tavg 1 
Rhow=FNRhow( Tavg > 

IF lstu=0 THEN 

Md = ( 6 . 7405 *Fn-H 3 . 0 2 ^ ) / 1 000. 

Md=Md* ( 1 .0365-1 . 96644E-3*T1 i +5 . 252E-6 * T 1 1*2 )/1 
ELSE 

Md= 1 . 04805E- 2 + 6 . 809 p 2E - 3 * F n 

Md=Md* ( 1 .0365-1 . 96644E-3*T1 i +5.252E-6*T1 r'2 >/ . 
END IF 
M f = Md/ Phow 
■J u = N f / ( PI ♦Dr 2/4 ; 

Vu5=Uw* ( D l / I . 27E- 2 ) 2 
Tcor=FNTfric< T2o ) 

T2o=Tco r 

Q=Md*Cpw* < T2o- T 1 i ) 

Op = Q/ < P I *Do *L ) 

I f t =0 

K uj = F N K w ( Tavg ) 

Muuj = FNMulj ( Tavg ) 

Rei=Rhow*Vw*Di/Nuuj * ASSUMED SAME FOR KORODENSE 
Pruj = FNPruj( Tavg ) 

Fe 1 =0 . 

Fe2=0 . 

Cf = 1 . 

Pruj f =Prw 
P e l f =Re i 
IF I n i =0 THEN 

Ome = Pei"Pe-.p*Pruj . 3333 *C f 
END IF 

IF I h i = 1 THEN 

Sra= . 88- ( . 24/(4. + Prwf ) ) 

Srb=. 333333^. 5*EXP( - . 6 * Prw f ) 

Qne= ( 5 . + .015*Reif*Sra*PrLjf‘'Srb ) 

END IF 

IF I h l =2 THEN 

Eps i = ( 1 . 82 *LGT ( Re i )- 1 . 64 ) ' ( -2 > 

P p 1 = 1 . +3 . 4 * Ep 5 l 

Ppk 2= 1 1 .7+1 . 8 * Prw '' ( -1 /3 ) 

Pp1=(Epsi/8 HRei*Pruj 

Pp2=( Ppk 1 +Ppk2*(Epsi/8 ) A . 5* ( Prw ” .6666- 1 ) ) 

Qroe = Pp 1 /Pp2 



.0037 

995434 
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247 9 END IE 
2 4 S I H i = f ud > D i • u i • u rn e 

2482 IF I - e = • J THEM 2 -jTO 2491 

2483 P I -P I * uo + bi 

2484 Ei 1 = 1 0 I - D l m P I ♦ 1 D i - 0 1 * . 5 

24 85 Ml = < H i * E 1 . •. k c ♦ H 1 1 ' . 5 

2435 P2= C I ♦ < Qi+G2 ' 

2487 A2 = < D2-Oi mPI + < Di+D 2 *.5 
2483 M2 = < Hi *P2 / ^ Kcu ♦ AZ ) > * .5 
2439 Pel =FNTanh( M 1 ♦ L 1 > / ( M 1 * L 1 ) 

2490 Fe2= s FNTann( M2*L2 )/(M2*L2 :■ 

2491 Qt=0/lPI*Di*(L^L1 *Fe1+L2*Fe2 * H i ) 

2492 IF I h i =0 THEN 

2 494 C f c = ; Muui/ FNMuw< Tavg'-Dt : )' . 1 4 

2495 IF AB5i< Cfc-Cf ;/Cfc / .001 THEN 

24 97 2 f = ( Cf +Cf - )* .5 

2500 GOTO 2461 

2501 END IF 

2503 END IF 

2504 IF I h i = 1 THEN 

2505 Pr'uj f c~FNPruj( Tavg + Ot ) 

250b Reifc=Ow*Di*FNRhow(Tavg+Dt >/FNMuw( Tavg+Ot ) 

2507 IF AES (( Prw f c -Pru f ) /Prw f c )>. 00 1 OR ABS( ( Re l f c-Re i f )/Re i f c ) . 00 1 THEN 

2508 Pruf = (Prujfc+Prwf)/2. 

2509 Reif=(Reifc^Reif)/2. 

2510 GOTO 2461 

251 1 ENO IF 
2513 END IF 

2516 Ift=Ifto 

2517 Lnt d=( T2o-T 1 i )/LOG( ( Tsteam-TI l )/lTsteam-T2o) ) 

2513 Uo=Q/ ( Lmt d *P I *Do *L ) 

2519 Ho* 1/(1 /Uo-Oo *L/ ( D l * * L + L 1 *Fe UL2 *Fe2 >*Hi ;-Rn ) 

2520 T cf= Op/Ho 

252 1 Cp5c=FNCpuj( ( Tconi-Tsteam)*.5 ) 

2522 H f g=FNH f g ( T 5 t earn ) 

2524 Tuio = Tsteam-Op/Ho 
2527 Tfilm=Tsteam/3+Two*2/3 
2530 Kf-FNKw(Tfilro) 

2533 Rho f =FNRhoui( T f 1 1 n ) 

2536 Muf=FNMuw(Tfilrn) 

2539 Hpq= . 65 1 *Kf*(Rhof’‘2*9.8l *Hfg/(Muf*Do*Qp ) ) A .3333 

254 1 Hnu3 5=.728*(Kf "3*9.81 *Hf g*Rho f s 2 / ( Muf *Do*Tc f > >*\25 

2542 A lp 1 = . 728*Ho/Hnuss 

2548 Tfm( J-1 ) = Tf llm 

255 1 Qpa( J- 1 )=Qp 

2554 Y=Hdq*Qp'' . 3333 
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2SS1 • =Qp 
2560 S.*. = S\ + < 

2562 d v = 5 / -*- y ' 



2566 Sv5=S,s+< 2 

256 5 5.*. y = z > +■ * * Y 

2572 01=500 

25” 01055=01/ ( 100-25 ,♦< Ts t ean-Troom ) l TO BE MODIFIED 

2575 Hf c-FNHf ( Tcon ; 

2534 Mdv=0 

2587 Bp = ( B vo 1 ♦ 100 .* ‘2/5.76 
2590 Hsc=Cpsc*(T5team“Tcon) 

2593 Ndvc= ( (Bp-Oloss /-Mdv*Hsc )/Hf g 

259b IF A B5< ( Mdv-Mdvc )/Mdvc ) .01 THEN 

2599 Mdv = (Mdv + ridvc,*.5 

2602 GOTO 2593 

2605 END IF 

2608 Mdv*- Mdv^Mdvc m .5 

26 11 Ug=FNUv s t < Ts t earn ) 

2614 Uv = Mdv *Ug/ A/. 

2 520 F = > 9.31 *Do *Mu f *H f g )/ ( Uv ’ x 2 * K f * ( Ts t earn- Two ) ) 

2623 Nu=Ho*Do/Kf 
2626 Ret-Uv*Fhof *Do/Muf 

2629 Nr=Nu/Ret A . 5 

2630 Hfuj = .96*( 9.31 *Hfg/Tcf ) * . 2 * K f . 8 * U v A . 1 *Rhof' s . 5/ ( Do *Mu f ) . 3 
2635 IF I ov = 2 THEN 

2645 PRINT 

2647 FRINT USING "SX," " Uu Rei Hi Uo Hfuj(DT) 



Hnui 0 ) 



2650 
,Hpq 

2651 PRINT 

2653 PRINT USING , 5X > “ “ Uv Ho q 

Hnu ( DT ) 

2655 PRINT USING " 5X , Z . DD , 1 X , 2 ( MZ . 3DE , 1 X ) , 2X , 3D . DO 




X ,3( MZ . 3DE , 1 X ) " ; Uv ,Ho ,0p 



Tcf 



NuRe 



F 



, Tc f , Nr , F , Hnu s s 
2656 PRINT 



2658 END IF 

2659 IF I ov = 1 THEN 

2660 IF I h i = 1 THEN 

2661 PRINT USING " 11 X ,DD , 2X , Z . DD , 1 X , 3< MO . 3DE , 1 X ) , 2 < 3D . DD , 1 X ) , D . DDD “ ; J ,Uw ,Uo , 
Ho f Qp , Tc f , Ts t earn , Sr a 

2662 ELSE 

2668 PRINT USING “ 11 X ,DD ,4X ,Z . DD ,2X ,2< MD. 3DE ,2X ) ,Z . 3DE ,3X ,30.00 ,2X ,3D . DD" ; J ,Uu» 
,Uo , Ho ,Qp , Tc f , Ts t earn 
2671 END IF 
2674 END IF 
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2575 


1 r _ m - 1 H £ hi 




:gb4 


i.* - hnd [j3d = i or iuii 0 then uutput 


OFiiep; Op , Ho 


2594 


F NO If 






^ — c 




I t_y 1 


w — 




2^08 


I- > = THEM 




2 2 5 


IF • I uj 1 i =0 HMD I j ob - 1 ) CP Ivjil = i THEN OUTPUT 


OF 1 lep ; Op ,Ho 


271 1 


INPUT 1 CHANGE TCCOL RISE 7 1 =Y f 2=N",Itr 




271 2 


If 1 1 r = ! THEN GOTO 2 384 




2713 


BEEP 




2 n 1 5 


INPUT 'OK TO STORE THIS DATA SET U=y i 0 = N)' 7, ‘ 


, 0 k 5 


2*71 5 


IF Q\ 5=1 T HEN 




2^25 


OUTPUT OF 1 1 e ; B v 0 I , B a n p ,Etp , F m , T 1 ,T2 , P v a p 1 


, P v a p 2 , E m f ( * ) 


2^35 


A i o 2 = A 1 p 1 ^H[p2 




2~4 9 


ELSE 




2752 


j=J- 1 




2"55 


END IP 




2^55 


BEEF 




2^5 1 


I NFU T ’Ul__ r hBP E 5E ANOTHER RUN 1 1 = / f 0=N;7 1 ' 


, G o _o n 


275 4 


Nrun= J 




2“ , 57 


IF 60 ..on 0 THEN Repeat 










^ / « 0 






-'7"7 


IF J Mrun THEN Repeat 




Z77S 


END IE 




Z779 


IF I j ob = 1 THEN 2809 




17Q7 


IF lwil=0 THEN 




2755 


h S 5 1 G N 0 F 1 1 e T 0 * 




2' 7 38 


I J 0D= 1 




2791 


I ujd= 1 




2797 


Im = 2 




2S00 


ASSIGN OF lie r 0 D_f 1 1 e$ 




2303 


GOTO 1 1 35 




2305 


END IF 




2309 


i- i f Q=a then 




23 1 2 


PRINT 




2 515 


91 = ' Nr ■ j r • S ■ > - 5 y * S < * ' ! Nr u n ♦ S >. s - S *. '2 




2 313 


Ac = ' 3y - 5 i * 5.* > / Nrun 




_ £ 2 *_ 


E^IMT USING "lOX, 1 " Least-Squares Line for Ho 


vs q curve." 


2524 


PRINT USING •’ 10X," " Slope = " " a M0.40E" 


5 SI 


2327 


PRINT USING " 1 0X , " " Intercept = “ ",MD.40E" 


; Ac 


2530 


END IF 




2333 


BEEP 




2343 * 


INPUT ENTER SANE TEMPORARY PLOT FILE NAME” ,Fp lots 


2353 


ASSIGN OFilep TO P_file$ 




2363 


FOR 1=1 TO Nrun 




2373 


ENTER OFilep;Qp,Ho 




2583 


Xc = LC'G( Op / Ho ) 
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2384 

:88s 

2885 

2337 

2838 

2383 

2390 

2391 

2892 

2893 
2394 
2895 
2998 

2897 

2898 
2399 
2902 
2905 
2908 
291 1 
2914 
2917 
2920 
2923 
2926 

2929 

2930 

2931 

2933 

2934 
2339 
2947 

2958 

2959 
2982 
2965 
2968 
2971 
2974 
2977 
2990 
2989 
2995 
299S 
3001 

A " ; J , 



rc = LOG( Op ) 

Z •. = Z • +- '< c 
Z • 2 = Z • 2 ♦ / c 2 
2 / = Z >. y + X c ♦ i c 
Z ,=Zy+X c 
.NEXT I 
Bb= . 75 

A a = E X P ( f , Z y - B b * Z /. ) / N r u n ) 

PRINT 

PRINT USING " 1 0X , ' ” Leas t -squares line for q = a*delta-T b 

PRINT USING "ISX.^'a = “ “ ,MZ . 4DE n ; Aa 
PRINT USING " 1 2X , " "b = , MZ . 4DE " ; Bb 

IF I f t = 0 THEN 
IF lpc = 0 THEN 
Qps=2.5E-5 

IF I i c - 0 THEN Hop = 9326 

IF Iic=l THEN Hop= 1 0 1 65 *(. 0 1 905/Do )'. 33333 
END IF 

IF Ipc=1 THEN 
Q p 5 = 7 . 5E + S 

IF 1 1 c =0 THEN Hop =7 1 76 

IF I lc* 1 THEN Hop = 7569*( .01 905/Do r.33333 
END IF 

Hos-Aa " ( 1 /Bb ) *Qp 5 ’( < B b — 1 ) / Bb ) 

IF lpc = 0 THEN A a s = 2 . 3 2 E J - 4 

IF I p c = 1 THEN Aas = 2.53E + 4 

A 1 p sm= . 876 • SWENSEN DATA 

IF lwil=0 THEN GOTO 2959 

Enrat=Alp2/Alpsm 

Enr =Ho5 /Hop 

PRINT 

FRINT USING "i2X ( n "Enhancement ratio (Del-T; 31 , 2D . 3D “ ; Enr a t 

ELSE 

PRINT 

IF I f t = I THEN 

Aas=2687.2 ' ZEBROLJSKI (U = 0.4 m/s; 

Aas = 2557.0*( .01 905/Oo )•' .33333 < UAN PETTEN (U = 0.25 m/s) 

END IF 

IF Ift-2 THEN Aas=9269.7*( .01905/Do) “.33333 
Edt =Aa/ Aas 
Eq = Edt ’’ ( 4/3 ) 

END IF 
BEEP 
PRINT 

FRINT USING " 1 0X , " *' NOTE : " " ,ZZ,"" data points were stored in file 

Q_f ile$ 
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3007 
3013 
301 B 
n " ; j 
3031 
3073 
3079 
3030 
3100 
31 16 
31 1 3 
3131 
3134 
3137 
3 1 30 
3133 
31 3B 
3139 
3143 
3 1 45 
3148 
3151 
3154 
3157 
31B0 
3163 
3166 
316 9 
3173 
3175 



. r 2 uere stored in data f.ie 



File 1 

i 1 g n 

'DUMMY" 
THEN 



TO 



CALL UilscmCi ) 



BEEP 
PRINT 

PRINT USING '10<, NOTE: X-v pa 

,PUt* 

EEEP 
rtSSIbN 
ASS I GN 
PURGE 
IF Iso= 

END 

DEF FNPvstC Tc > 

CDM /Fid/ I ft , Istu 
DIM K ( 3 ) 

IF I f t =0 THEN 

DATA -7.691 23456 4 ,-25.08023696 ,-163.1 706546 ,64.23235504 ,-113.9646225 
DATA 4.1671 I 732 ,20.3750676 ,1 .E9,6 
READ K< * ) 

T=( Tc + 273 . 15 ■ / 647 . 3 

Sun = 0 

FDR N=0 TO 4 

Sum=Sum+K ( N ) ♦ ( 1 -7 ) " ( N+ 1 ) 

NEXT N 

Br=Sum/< T* < I +K( 5 ) * ( 1 -T ) + K( 6 ) ♦ U -T ) "2 ) )-< 1 -T )/ ( K ( 7 )» ( 1 -T ) A 2+K ( 8 ) ) 
Pr=EXP( 6r ) 

P = 22 1 2 0 0 0 0 * P r 
END IF 

IF I f t = 1 THEN 

T f = Tc* I .8 + 32 + 459.6 

P=10 / '( 33.0655-4330. 93/ Tf -9. 2635 »LGT< Tf ; + 2 . 0539E-3»Tf ) 

P=P* 1 0 1 325/ 14.696 



’0 



3178 END IF 

3131 IF I f t =2 THEN 

3134 A = 9. 394685-3066. 1 / ( Tc+273 . 1 5 ) 

3187 P= 1 33 . 32 * I 0 A 

3190 END IF 

3193 RETURN P 

3196 FNEND 

3199 DEF FNH f g( T ) 

3202 COM /Fid/ Ift .Istu 
3205 IF I f t =0 THEN 

3203 Hfg=2477200-2450*(T-10) 

321 I END IF 

3214 IF Ift-1 THEN 
3217 Tf=T*l .8+32 

3220 Hfg=7.0S57857E+l-Tf*(4.83S052E-2 + l . 26 1 9048E-4* T f ) 

3223 Hfg=Hf g»2326 
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3226 

3229 

32 j2 

3235 

3236 

324 I 

3244 

3247 

3250 

3253 

3256 

3259 

3262 

3265 

3263 

3271 

3274 

3^77 

3230 

3233 

3286 

3289 

3292 

3295 

3298 

3301 

3304 

3307 

3310 

33 1 3 

3316 

3319 

3325 

3323 

3331 

3334 

3337 

3340 

3343 

3346 

3349 

3352 

3355 

3358 

3361 



END IF 

IF I f t =2 THEN 
74 =T+273. 15 

H f g = 1 . 35264E + 6 - T 1 * ' 6 . 38263E+2+Tk * .T’47462 
END IF 
RETURN Hfg 
FNEND 

DEF FNMuw( T > 

COM /Fid/ I ft , Istu 
IF I f t =0 THEN 

A = 247 . 8/ ( T+ 1 33 . 15 ) 

Mu = 2 . 4E-5 * 1 0"A 
END IF 

IF Ift=l THEN 

Mu =8. 962981 9E-4-T * ( 1 . 1 0 94609E-5-T *5 . 56682 9E-9 ) 

END IF 

IF Ift«2 THEN 

Tk = l / ( T + 2" 1 3. 15 ) 

Mu=EXP ( - 1 1 .01 79 + Tk ♦ ( 1 .744E + 3-U ♦ ( 2 . S0335E + 5-Tk *1.1 26B 1 E + 8 ) ) ) 
END IF 
RETURN Mu 
FNEND 

DEF FNUvs t < Tt ) 

COM /Fid/ Ift , I s t u 
IF Ift =0 THEN 
P=FNPv 5 t < Tt ) 

T = T t +273 . 15 
X= I 500/T 

FI =1 /( 1 + T * 1 .E-4 ) 

F2 = ( 1 -EXP ( -X ) ) ' 2 . 5*EXP< X )/X“ .5 
B=.001 5 *F1 - .000942 *FZ- .0004882 *X 
K = 2 * P / (461 . 52 ♦ T ) 

U = ( 1 + ( 1 ^ 2 ♦ B * K )•' .5 )/K 
END IF 

IF I f t = 1 THEN 
T f = T t * 1 .8 + 32 

U=1 3. 955357- Tf *( . 1 6 1 27262- T f *5 . 1 726 1 90E-4 ) 

U=U/ 16.018 
END IF 

IF I f t =2 THEN 
Tk =Tt +273 . 1 5 
P = FNPvs t ( Tt ) 

U= 1 33 . 95 * Tk /P 
END IF 
RETURN U 
FNEND' 
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536*1 EF ~ ' 1 C 3 ui 
3367 COM F : a ' I •' t , I 5 1 u 
i F I f t =*0 T HEN 

3 37’3 Co -j = 4 .3 ! 1 3 Z 353 -T. ... 366 366-3-1 •' 4 . 4 3 3616-5^3 . 71 428E-7*T 

3375 £MD IF 

3379 ’ p 1 ►*. = 1 — h£ N 

3363 Cpu = 9 . - 5073 3E-1/'T*(9.o4004o - E-4'*'i,7207~ , 3 — t-b* 

3365 END IF 

3333 IF 1 1 1 = 3 THEN 

3391 TF = T + 373 . 1 5 

3394 Cpw = 4. ’ S 6 S * ■ 1 . S834E-3 + TI' »< I. 350S3E-3-TI- m 7 . 224E-G-U *7 . 6 1 748E-9 > > 

3397 END IF 

340(3 PETIJRN Cpw*'0O0 

3403 FNEND 

3406 DEF FNFhCuj' T 1 

3439 COM , Fid/ If t , istu 

3413 IF i f t =0 THEN 

3415 Ro = 999.53 946 + T*< ,0i 36 3-T*' 5. 4335! 3E-3-T* 1 .334U7E-5 ) ; 

3418 END IF 

3431 IF I f t = 1 THEN 

34 34 Ro =1 . 5307479E + 3-T » ' 2 . 3 1 86346 +T * 3 . 357339 1 E-3 ) 

3427 END IF 

3430 IF I f t = 2 THEN 

3433 T1 =T+273 . 1 5-338 . 1 5 

3436 Df = 9.24848E-4fTkn S . 27 96E-7 + TT *• < 9 . 2 444E- 1 0 + Tk »3.057E-1 2 ) ) 

3439 F, o = I / U f 

3442 END IF 
3445 RETURN Ro 
3448 FNEND 
3451 DEF FNPru(T) 

3454 Pr 'j=FNCpu ( T ) »FNMuw< T )/FNKw( T ) 

3457 RETURN Fry 
3460 FNEND 
3463 DEF FNKw(T) 

3466 CON /Fid/ I ft , Istu 

3469 IF I f t = 0 THEN 

3472 / = ( Ti-273. 15 )/273. 15 

3475 Kw = - . 92247 + X » < 2 . 8395-X*( I . 8007-X*( . 52577- . 07344 * X ) ) ) 

3478 END IF 

3481 IF I f t = 1 THEN 

3484 Ku=8.2095238E-2-T»‘- 2 . 22 1 4286E-4 + T *2 . 3809524E-8 ) 

3487 END IF 

3490 IF I f t = 2 THEN 

3493 Tk-T+273.15 

3496 Kui=4. l868E-4*(519.442+.320920*Tk ) 

3499 END IF 
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3502 RETURN Kw 
3505 FNEND 
3538 OEF FNTanhi x 
35 I I P=E*P( X ) 

3514 Q=E <P ( ' 

3517 Tann=( P-Q )/( P + Q ) 

3530 RETURN Tanh 
3523 FNEND 
3526 DEF FNTvsvOJi 
3529 CON /Cc/ C(7 ) 

3532 T=C ( 0 ) 

3535 FOR 1=1 TO 7 
3538 T = T + C( I >*()” I 

3541 NEXT I 
3547 RETURN T 
3550 FNEND 
3553 DEF FNHflT, 

3556 COM /Fid/ I ft ,Istu 
3559 IF lft=0 THEN 

3562 Hf = T»( 4. 203849 -Tm 5.881 32E-4-TM .551 6031 7E-6 ) ) 

3565 END IF 

3568 IF I f t = ! THEN 

3571 Tf =T» I .8+32 

3574 Hf =8 . 207557 1 +T f ♦ ( . 1 9467857 + Tf * 1 .32 1 4 286E-4 ) 

357" Hf =Hf *2 . 325 

3590 END IF 

3583 IF I f t =2 THEN 

3586 Hf =250 ' TO BE UERIFIED 

3589 END IF 

3592 RETURN Hf • 1000 

3555 FNEND 

3598 DEF FNGradl T ) 

3601 5rad=37 . 9853 +. 1 04388 ♦ T 

3604 RETURN Grad 
3607 FNEND 
3610 DEF FNTvsp ( P ) 

3613 Tu= 1 90 
3616 T 1 = l 0 
3619 Ta = ( Tu + TL )* .5 
3622 Pc = FNPvst( Ta ) 

3625 IF ABS( (P-Pc )/P .0001 THEN 

3628 IF Pc<P THEN Ti=Ta 

3631 IF Pc>P THEN Tu=Ta 

3634 GOTO 36 1 9 
3637 END IF 
3640 RETURN Ta 
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3643 


FrjEND 




bo 4 b 


QEF P N8 i groa ■ T > 




bo 4 3 


<=S4^. 3- r -.:^3 . ; 5 




665Z 


3 - . * 1609^6^07' 1 1 ^ . 3 3 * ^ 1 + i . 12l404688E-3~5 


. 752305 1 8E-5 • x + 1 . 


. 1 4 97 


1 523E-I i m 3 




6655 


RETURN £*.001 * a 2 




6659 


FNENO 




6326 


SUB to l i 5 on ( Ci > 




6 31 9 


COM Ul i 1 / Nrun.It.'b.Iwth^mc.Ife.Ijob.Iud 


, I fg , Ipco , I f to , I ui 


zu ,Re 


A p , R m , A a 




6332 


COM /Fid/ I f t , I s t u 




6S33 


COM . 6eom/ D1,D2,Di,Do,L,Ll f L2 




6334 


COM /Eric/ I 5 t uo , Inn , Iu t ho ,Uu ,Uws , I narno 




6836 


DIM Emf ( 20 ) ,Bvo< 42 ) ,Bam( 42 ) ,Eata( 42 ) ,Ear( 42 ,7 ) f Fma( 42 ) ,T 1 


6345 


IF I oc “0 THEN 




6847 


PRINT USING " G/ , ' " Nu 5 se 1 1 theory is 


used for Ho 


6348 


ELSE 




6349 


PRINT USING ' 1 6X , " Fu j 1 1 correlation 


used for Ho 


6350 


END IF 




6853 


BEEF 




6856 


INPUT "RE-ENTER DATA FILE BEING PROCESSED" ,D_f 1 le* 


6859 


BEEP 




6862 


INPUT “GIUE A NAME FOR XY PLOT-DATA FILE 


" .Plot* 


6365 


CREATE BOAT Plot* ,10 




6368 


ASSIGN @Io_path TO Plots 




687 1 


J j=0 




6874 


ASSIGN AFile TO D_file* 




6877 


ENTER ®Fi 1 e : I f g , Inn 




6 o 7 8 


IF I s t u - 0 THEN 




6383 


ENTER @Fi le ; Ddd ,Ddd ,Ddd ,Ddd 




6384 


ELSE 




6585 


IF I f g = 0 THEN ENTER @F i 1 e ; I ui t 




6 386 


IF 1^=1 THEN ENTER @F l le i Fp ,Fw ,Fh 




6387 


END IF 




6388 


IF Jj=0 THEN 




6395 


IF lhi=0 THEN C i = . 027 




6896 


IF IHii = 1 THEN C l = 1 .00 




6897 


IF I h i =2 THEN C l = 1 .00 




6899 


IF I f g=0 THEN Alp=1 .2 




6900 


IF Ifg=l THEN Alp=2.S 




6901 


IF I f t =2 AND I f g= 1 THEN Alp=5.0 




6904 


END IF 




6907 


J=0 




6910 


S A = 0 




6913 


Sy = 0 




6916 


S>:5 = 0 





113 



6319 : A / =0 

6922 ' PEnD DATA FROM A USER-SPECIFIED FILE IF INPUT MODE In; = 2 

5 925 IF j j — 0 Then 

6 326 IF I s t ■ j - 0 ThEN 

6 931 ENTER ^Fiie;6voi,0anp,Etp,Fn,Tl,T2 > Ddd,Ddd,Enf<*' 

5932 ELBE 

6934 ENTER , ?File;6voi l Ganp > Utran,Etp,EHf(0),Emf(i) l EHf(2) J Emf<3) J Emf(4 f F 

n , T I , T2, Phg, P water 
6936 END IF 

6933 5 vo ( J ; = Bvo I 

6939 Eam J ;=5anp 

6940 Eata(J;=Etp 

6943 Ear ( J ,0 >=6nf ( 0 ) 

6946 Ear ( J , 1 > = En f ( 1 ) 

6949 Ear>.J,2.> = Enf(2) 

6952 Ear' J , 3 > = En f ( 3 ) 

6955 E ar ( J , 4 ; = E n f ( 4 ) 

6956 IF I 5 1 ii ~ 1 THEN GOTO 696 1 

6958 Ear< J ,5 )=Enf < 5 ) 

6959 Ear(J l 5)=Emf(6) 

5960 Ear( J ,7 )=Enf ( 7 ) 

696 1 Fna( J ) = Fn 

6962 T 1 a ( J > = T 1 

6964 T2a ( J ) = T2 
6967 ELSE 

6 970 8vo 1 =B vo( J ) 

6 373 6ann=6an( J ) 

6976 E t p = Ea t a ( J ) 

6979 Enf ( 0 )=Ear' J ,0 ) 

6982 Enf ( 1 ) = Ear< J J ) 

6965 Enf ( 2 )=Ear< J ,2 5 

6988 Enf ( 3 )=Eam J ,3 ) 

6 991 Enf ( 4 ) = Ear ( J , 4 ) 

6992 IF I s t u= I THEN GOTO 6997 

6994 Enf ( 5 )=Ear( J ,5) 

6395 Enfv6) = Ear(J,6) 

6996 Enf ( 7 )=Ear( J ,7 ) 

6997 Fn=Fna( J ) 

6998 T 1 -T 1 a ( J ) 

7000 T2=T2a ( J ) 

7003 END IF 

7004 IF I s t u = 0 THEN 

7006 T sat=FNTvsv57( < Enf < 0 ) + En f ( I ) >/2 . ) 

7007 Taat =T sat-273 . 15 

7009 Ti=FNTv5v58( Enf ( 3 ) ) 

7010 To1=FNTvsv58(Enf (4 ) ) 
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”3 I 3 T:=t l -^73. ' 5 

"013 r 3i-T 0 !-:-3. !5 

. 0 1 5 ti f p I =Em^ 3 f E * p 30. 

”015 -ittie-E.ESd’L-J - 1.5246 4 E-b*Etp +- 1 . _ \ 7 0 ] E 0 * b t p 1 3~5.ilt4t:-15*E t . pi 

3 0 I E - 1 i 3 * b t p ! 4 

7017 7r : 5=0 r . pe ♦ E : p / 10. 

7 0 ! 3 7 c = 7 i+Tr i 5 

7019 ELSE 

7030 7sat=FN7ysv< E ^ f •' 0 m 

”031 7 i = F N T j s v v E f ' 3 ) ) 

7033 Grad=FNGrad( (71+73 « . E j 

7033 7o = T i +mBS< E t p > / M O^Graa - ♦ . E + b 

”334 END IF 

7035 ' 7ELCULE7E "HE LOG-HErN- T E M PE p H TI JRE DIFFERENCE 

”036 I c I 5 t u = 0 H N D 1 1 t \ = 0 7 n E N 

”037 r - = 7 i 

”338 T 1 = 7 c 1 

”339 END ;F 

7030 IF 1 1 m = 1 THEN 

7031 7 f =-7 1 

7033 71=73 

7033 END IF 

7034 IF 1 1 m = 3 7 HEN 

7035 7 f = 7 i 

7036 71=13 
7038 END IF 

”039 7 a v g = ( 7 f <- 7 1 *.5 

7040 Tr L 5e = 7 1 -7 f 

”041 L' y itd = 7rise/L0G ( v 7 s a t - 7 c : 7 s a t - 7 1 )) 

7043 I f t = 0 

”043 Dpuj = FNCpw( Tavg > 

*^044 Rhaui = FNRhcw'7a/g) 

”04 5 H % = FNKuj ( 7a > 

“ 7 048 Nuoja = FNMuw( 7a vg ) 

”05 1 p r =F NP^ w ( 7 a v g ; 

”054 Ift-Ifto 
7055 IF I 5 t u =0 I HEN 

”05” Mdt = ( 6 .7409*Fm+1 3.037 > / | 0 0 0 , 

”OE0 Md=Nd t * ( 1 .0365-7 f*< 1 . 96 6 44E-3- 7 f *5 . 3S3E-6 ) )/\ .0037 
”061 ELSE 

7063 Mdt = 1 .04S05E-3f6 . 30933E-3 *F-i 

706 3 Md=Mdt * ( 1 .0365- 7 f *( 1 . 96644E-3-7 f *5 . 353E-5 ) >/ .995434 

”065 END IF 

7066 V f =Md/Rhou 

7067 0uj = Vf/< PI *0i ■ 3/4 ) 

”065 Jws^wmDi/1 . 37E-3 ) 3 
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^06 9 Tcor=FNT f r i c ' Tr 1 se * 

7070 rri5e=Tccr 
7 0 7 1 T c o n = F N T v 5 v 5 8 ' E m ? ( 7 

"07E Tcon=Tcon-273. 15 
7074 Cpsc = FNCpuu f rcnrT^at *.5.' 

7076 Hfg=FNHfg / T 5 a i > 

7073 01=500 

7030 Qloss=Ql/U00-25 >*< Tsat-Troon > « TO 5E MODIFIED 

7082 M d v = 0 

7084 5p = < Bvo 1 * I 00 > ,s 2 /5 . 76 

7085 Hsc = Cp-5G*( Tsat-Tcon ) 

7036 Mdvc = ( ( Bp-Q 1 os s )-Mdv*Hsc H f g 

'7086 IF ABS ^ (Mdv-Mdvc ) / Mdvc ) - .01 THEN 

7050 M d v = ( Mdv+Mdvc )• .5 

7052 GOTO 7096 

7094 END IF 

7096 Md*=< Mdv+Mdvc . 5 

7098 Vg=FNVvs 1 ( Tsa t ) 

7100 0 v = M d v ♦ 0 g / A x 

7108 Q=Md*Cpw*Trise 
71 1 1 Qp = Q/ (PI»Do*L ) 

7114 Uo=Qp/Lmt d 
7117 Re=Rhow * D i / M uwa 
7120 Fe 1 =0 
7123 Fe2=0 

7126 C f = 1 . 

7127 Prwf = Prui 

7128 Reif-Re 

7129 lft=0 

7131 Two = Tsa t -5 

7132 Tfiln=T5at/3+Two*2/3 
7135 h f = F N K oj ( T film ) 

7138 R ho f =FNR hou T f 1 1 m ) 

7141 Muf=FNMuw(Tfiirn) 

7144 Hfgp = FNHfg( Tsat >+.6S*FNCpu( Tfilm;*(Tsat“Two ) 

7147 1 New = K f * i Rho f '* 2 * 9 . S 1 *Hfgp/(Muf*Do*Qp ) ) ' .3333 

7148 New= ( *■ f “ 3 * 9 . 8 1 *H f gp *Rho f *' 2 / ( Mu f *0o * K T sa t -Two ) ) ) .25 

7150 IF I oc = 1 THEN 

7153' New=K f * ( ( 9 . 8 1 *H f gp / Qp ) * . 25 ) ♦ ( ( Mu f *0o ) * < - . 375 ) >♦< Rhof * .625 )*( Uv .125) 
7154 Ne w = ( 9. 81 *Hf gp / ( Tsat-Two ) ) A . 2 *K f A . 8 *Uv *’ . 1*Rhof A .5/(Do*Muf ) A .3 

7156 END IF 
7159 Ho=Aip*New 
7162 Twoc=T sat -Qp /Ho 

7165 IF ABS ( ( Twoc- Two )/Twoc ) ' . 00 1 THEN 
7168 Tuo=Tuoc 

7171 GOTO 71 32 
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7174 END IF 

"'175 ^e«Di=Fevp 

^"3 4 I F I h i = 0 then 

7185 8^ega = Fe R e • p i ♦ p y .3353* f 

7 1 3 T END IF 

7183 IF I h i = 1 'HEN 

7189 S r a = . 3 3 - ( .24/' 4. -•-Pry f ); 

7190 SrD=.333333+.5*EXPC-6 # Prwf ) 

7191 Qnega=(5.+.015»Rei* Sra*Prwf 5 r P ; 

7193 END IF 

7193 IF I hi =2 THEN 

7194 Ep 5 l = ( 1 . 9 2 ♦ L G T ( Re j - 1 . 54 > •' -3 ) 

7195 FpH s 1.+3.4*Ep5i 

7196 P p k 2 = 1 1 .7+1 . 8 ♦ Pry ' ( -1 3 ) 

7197 Fp1=(Ep5i / 3/*Pe*P^'j 

7198 Pp2 = < Ppl- 1 + Pp1 2 ♦ < Epsi '3 ; .5* Pry . 6 6 6 6 - 1 ) ; 

7199 Ornega = Pp 1 / Pp2 

7200 END IF 

7202 Hi=Kw/Di *Ci ♦Onega 

7203 IF I f e = 0 THEN 72 1 5 

7204 P 1 =PI ♦< Di +D1 ) 

7205 P2 = P I * ( Di +D2 ) 

7206 A1=(D1-Di .♦PlMDi^DI *.5 

7207 A2 = ^ D2-D i ) *PI ♦ < Di +D2 ) * . 5 

7208 M 1 - ( H i ♦ P 1 /< Kcu*nl i > ' .5 

7209 N2 = ( Hi * P 2 / 1 ^ c u ♦ H 2 ) ) ' .5 

7210 Fei=FNTanh^M1*LlV(N!*Ll) 

7213 Fe2=FNTanh( M2 *L2 )/ v M2 *L2 ) 

7216 Dt=Q/\PI*Di*(L+Ll *Fe1 + L2 *Fe2 ) *H i i 

7217 IF lhi=0 THEN 

7219 Muuj i =FNMuw ( Ta vg + D t ) 

7222 Cfc = ( Muwa/Muwi ) ' . U 

7225 IF nS3UCPc-Cf »/Cfc ^ .001 THEN 

7223 Cf=iCf+Cfc)*.5 

7231 GOTO 7184 

7232 END IF 

7234 END IF 

7235 IF I h l = 1 THEN 

7236 Prwfc=FNPry( Tavg+Dt ) 

7237 Pei fc=Ou*Di»FNRhow' Tavg + Dt )/FNMuw( Tavg + Dt ) 

7239 IF ABS( (Pryfc-Pruf //Pruifc >• .001 OR ABS( (Peifc-Reif >/Rei fc >/ .001 THEN 

7240 Prwf =( Prwfc+Prwf >/2. 

7241 Reif=(Reifc+Reif)/2. 

7242 GOTO 7184 

7243 END IF 
7245 END IF 
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7246 Ift=ifto 

7247 k =0o * New * L / ( Gmega*Ku*'- L + L ! * F e 1 + L 2 * F e 2 ) ) 

7243 > = N e w # ( 1 / U o - p m ) 

7249 1 COMPUTE COEFFICIENTS FOR THE LEAS T- SQUARE S-F I T STRAIGHT LINE 
'7250 IF Jp=J THEN OUT p, JT @Io_pa t h ; X , Y 
7252 5,.=S^< 

72E5 S/=Sy+r 
~ 1 5 8 S/s = S,*'S + X*X 
"'261 S y = S a y + X * Y 

7264 IF Im=l AND Jj=0 THEN OUTPUT <§File;6vol J Barnp,Etp,Frn,Tl } T2,Pvapl,Rvap2,Em 

* ) 

7267 J = J + 1 

7270 IF J Nrun THEN 6925 

7273 Sl = (Nrun + Sxy-Sy*S.* ) / ( N r u n * 5 * s - S * * 2 ) 

7276 IF I w i 1 = 2 THEN 
"~7 IF 1 h i = 0 THEN 51 = 1/ .027 

9 IF I h i = 1 THEN S I = 1 '1.00 

7289 IF I h l = 2 THEN SI = 1 / 1 .00 

7291 END IF 

72 94 Ac = ( Sy-S 1 *S < ) /Nrun 

7297 Cic=l/Sl 

7 300 Alpc=1/Ac 

7303 Jj=Jj+I 

7306 IF Jp=1 THEN Jp = 2 

7309 Cerr = A6S ( f C l c -C l )/C l c ; 

7312 Aerr=ABS ( < A i pc -A 1 p ; / A 1 pc ) 

7315 IF Cerr .001 OR Aerr .001 THEN 

7313 C l = ( C i c +C l ) ♦ . 5 

7321 A 1 p = ( Alpc+Alp )* .5 

7324 BEEP 

7327 IF I j ob- 1 THEN 6907 

7330 ELSE 

7333 IF J p =0 THEN Jp=1 

7336 END IF 
7339 IF Jp=1 THEN 6374 
7342 Ci=<Ci+Cic m . 5 
7345 PRINT 



7346 


IF lhi=0 


THEN 












7348 


PRINT 


USING 


" 1 0X , 


"Cl 


( baaed on 


Sieder-Tate ) = 


,Z.4D" ;Ci 


7349 


END IF 














7350 


IF I h i = 1 


THEN 












735 1 


PRINT 


USING 


" 1 0X , 




( based on 


S 1 e l c her-Rouse ) 


,Z.4D" 5 Ci 


7352 i 


PRINT 


USING 


" 10X 4 


, " "Re 


exponent 


for S leicher-Rouse = 


, D . ODD " ; Sra 


7353 


END IF 














7354 


IF I h l =2 


THEN 












7355 


PRINT 


USING 


’ 10X , 


"Ci 


(based on 


Pe t u\r hov-Popo v ) = 


"" ,Z.4D" iCi 
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"356 


nNu I c 










"357 


IF I dc -0 r u c rxj 










7 353 


P=INT j 5 I N-j 0" 


, Hip na based cn Nusse 


it ' Tde 1 1 = , 2 . 4Q " ; A lp 


^353 


ENG IF 










"36 3 


IF Icc=< 'HEN 










"331 


F-IN T 'J = ING 0X 


, ’ A I p n a based cn huj i l 


[ T d e i ' ) = , 2 . 4 D " ; n l p 


"-63 


ENG IF 










"363 


IF I n a m = 5 jR I n a m = 6 


THEN 








"364 


IF I h i = 0 THEN 










^365 


IF I p zo -0 AND 


I n n = 0 


THEN 


Alpsn= .8318 


'NO INSERT , VACUUM, S-T 


"367 


IF Ipcc=1 AND 


I nn = 0 


THEN 


A 1 p sm= . 7793 


i NO INSERT .ATMOSPHERIC ,S- 


< j 6 2 


IF I pc o = 0 AND 


I nn = 3 


THEN 


A lp sm= . 7854 


'HEATEX .UACUUM ,S-T 


7369 


IF Ipco=1 AND 


I nn = 3 


T HEN 


A i p sn= . 7769 


'HEATEX .ATMOSPHERIC ,5-T 


7371 


END IF 










7 g7 *_ 


IF I n i = 1 i nEN 










7373 


IF I p c o = 0 H N D 


I n n = 0 


THEN 


A 1 p sm= . 66 1 3 


'NO INSERT .VACUUM ,S-R 


"3"4 


IF I p c o = i A N D 


I Hn = 0 


THEN 


A i □ 5H= . 83 1 8 


'NO INSERT .ATMOSPHERIC ,S- 


"375 


I c I p c c - 0 hND 


I nn = 3 


THEN 


Aipsm=.7791 


' HEATEX .VACUUM ,S-R 


"376 


IF I pco= 1 AND 


I nn = 3 


THEN 


Aip 5 m=.7939 


'HEATEX .ATMOSPHERIC ,S-P 


73" 3 


END IF 










7379 


IF Ihi=3 THEN 










7360 


IF Ipco-0 A ND 


I nn = 0 


THEN 


A 1 p sm= . 8305 


'NO INSERT .VACUUM , P-P 


7361 


IF I p zo= 1 AND 


I nn = 0 


THEN 


A 1 p am= . 7654 


'NO INSERT .ATMOSPHERIC ,P- 


7383 


IF lpco=0 AND 


I nn = 3 


THEN 


A 1 p 5n= . 7670 


' HEATEX .VACUUM ,P-P 


"333 


IF I p c o = 1 mND 


I nn = 3 


THEN 


A 1 p sm= . "708 


'HEATEX .ATMOSPHERIC .P-P 


7385 


END IF 










"366 


END IF 










7387 


IF Inam=4 THEN 










"390 


IF I p c c = 1 THEN A 


1 p srn= . 


876 


• SWENSEN DATA 


BASED ON DEL-T 


"391 


END IF 










"393 


IF lnam = 0 OR I n a m = 3 


THEN 








"393 


IF lpco = 0 THEN A 1 p s m = . 


53 


MJP Ml S T 0 1 0 3 




7334 


I f I p co= 1 THEN Alp5H=. 


68 


5MTHSTA65 




7337 


IF Ift = 1 THEN h i p 5n= . 6 


77 


i VAN PETTEN ( V = 0.25 n/s ) 


7399 


IF I f t =3 THEN Alp5n=1. 


262 






"399 


END IF 










740 1 


IF I nam= 1 THEN 'MITROU ALPHA 1 


FOR P-P FROM 


REPROCESSING 


"403 


IF lpco = 0 THEN A 1 p s m- . 


8437 






7 403 


IF Ipco=1 THEN Alpsm=. 


8418 






7404 


END IF 










7405 


Et=Alp/Alpsm 










7406 


Eq=Et 1 .333333 










7407 


PRINT USING " 1 OX ,""l 


Enhanc emen t 


(q ) 


= ■* .0D.3D" ; Eq 


7408 


PRINT USING “ 10X , " "1 


Enhancemen t 


( Del-T ) 


= ,DD . 3D " sEt 


7409 


ASSIGN @File TO * 










7410 


SUBEND 
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7721 DEF FNTv5v55( V ; 

^731 COM /Cc55 / T55( 5 ) 

7"U1 r=T55(0) 

"751 ROR 1=1 TO 5 

7761 7 = T + T 5 5 ’ I ;*U' I 

7771 NEXT I 

7781 RETURN T 

7791 FNEND 

7801 DEF FNTv 5 v56 ( U ) 

781 1 COM /Cc56 / T56( 5 ) 

7321 T = T5 6(0) 

7831 FOR 1=1 TO 5 

784 1 T = T + T56 ( I )*V I 

7351 NEXT I 

7861 RETURN T 

7871 FNEND 

788 1 DEF FNT , 5 v57 ( U t 

7891 COM / Cc57 / 157 ^ 5) 

7901 T = T57 ( 0 ) 

791 1 FOR 1=1 TO 5 

7921 T = T + T57 ( I )*U ,S I 

7931 NEXT I 

7941 RETURN T 

7951 FNEND 

"961 DEF FNT.5v58<U) 

7971 COM / C c 5 3 / T53;5 ) 

7981 T=T58 (0) 

7991 FOR 1=1 TO 5 

8001 T = T + T53( I )*V I 

8011 NEXT I 
8021 RETURN T 
8031 FNEND 

8100 DEF FNTf~icl Tcor ) 

8110 COM /Fric-' I 5 tuo , I nn , Iut ho ,Uuj 5 , 1 narno 

8115 IF I s t u o = 0 ANO I lj t ho=0 THEN 1 SWENSEN FRICTION FAC. FOR COPPER TUBE 

8120 IF I nn = 0 AND Uui^.S THEN Tcor = Tcor-< -2 . 73E-4+ 1 . 75E-4*Uw + 9 . 35E-4*Uu "2- 1 . 9 

6E-5 ♦Uu " 3 ) 

8130 IF I nn= 1 THEN Tcor=Tcor-< -6 . 44E-5+ 1 . 7 1 E-3 *Uuj + 4 . 45E-4*Uuj ' 2 + 4 . 07E-5*Uu '3 ) 

8140 IF I nn=2 THEN Tcor = Tcor-( -3 . 99E-4+2 . 75E-3*Uu+l . 45E-3*Uu *2 + 8 . 1 6E-5 *Uuj a 3 ) 

8150 IF I nn = 3 THEN Tcor=Tcor- ( 8 . 57E-5+1 . 23E-3*Vw+ 1 . 08E-3*Uw A 2+8 . 1 GE-5*Vw*3 ) 

8160 END IF 

8170 IF I s t u o = 0 AND Iutho=2 THEN 1 OKEEFE FRIC. FAC. FOR SMOOTH TITANIUM TUBE. 
8180 IF I nn=0 ANO Uu>.5 THEN Tcor=Tcor-< -4 . 62E-5-7 . S3E-4 *Uuj + 1 . 80E-3*UuT2-8 . 8 

4E-5 * Uw ' 3 ) 

8190 IF Inn = 3 ANO Inamo=5 THEN Tcor = Tcor-( 2 . 09E-4+9 . 74E-4*Uw+2 . 1 2E-3 «Uuj A 2~3 . 

31E-5*Uui A 3 ) 
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lu cn u) 



3 1 ! 3 1 IF Inn-3 AND Inano=6 THEN Tcor = Tccr-< 1 . 9555E-4 + 3 . 972 1 E-3 *Ulj + 3 . 1 27E-4 * Uw 

2+3.5! 9E-4*Uw 3 > 

200 END IF 

213 IF I s t u j ~ 0 AND I *jj t m o = 1 THEN ‘FRICTION FACTORS FOP LFD KOFODENSE 
310 I p I^n-O AND Jui .5 THEN Tcor= Icor-( -3 . 386E-4 + 1 . 88E-3 *Uu+6 . 0 I 3E-4 3-4 

. } 3IE-5 ♦ ,“uj 3 ' 

6330 IF Inn = 3 THEN Tcor=Tccr-( 3 .089E-4+9 . 202E-4 *Uw+ l . 893E-3 " 2 -2 . 78 1 E-5*vw 

'3 ) 

8240 END I F 

8342 IF 1 5 1 u 0—0 AND Iwtho-5 THEN ‘FRICTION FACTORS FOR NHT KORODENSE 

8244 IF I nn=0 AND Uuj ; . 5 THEN Tcor = Tcor- ( E . 7 9E-5+ 1 . 632 E~3 * Vw+8 . 409E-4 * Uw ” 2+1 . 

1 1 1 E-4 +i)v 3 ) 

8245 IF Inn = 3 T HEN Tcor=Tccr~( 2 . 564E-4+6 . 263E _ 4*Uui + 2 . S03E-3 ,N 2+7 . 830E-6*Uw 
3 ) 

8247 END IF 

8250 IF Istuo-i THEN 

8250’ IF lnn = 0 THEN Tccr=Tcor-( . 0 1 38+ . 00 1 *Uw‘'2 ) 

8270 IF I nr= 1 THEN fcor- Tcor- . 004 *Uws ' 2 

8280 IF Inn = 2 THEN Tcor-Tcor- . 004 + Uws ,N 2 

S2S0 END IF 
8300 RETURN Tcor 
8310 FNEND 
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